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ABSTRACT 


a  computer  program  has  been  developed  for  one-dimensional  nonequi librium 
reacting  gas  flow.  The  program  is  written  in  Fortran  IV  and  is  compatible 
with  the  IBM  7044/70%  direct  coupled  di gitai  computer  system  at  Wright- 
Fatterson  Air  Force  Base,  Ohio.  In  addition  to  nonequilibrium  chemistry,  the 
program  includes  nonequilibrium  'idbrational  and  electronic  energy  relaxation 
and  coupling  effects  between  these  energy  modes  and  xhe  chemistry.  The  form¬ 
ulation  is  based  on  a  one-dimensional  flow  matching  either  a  prescribed  pres¬ 
sure  or  area  variation  along,  a  streamtube.  Thermodynamic  properties  are 
computed  by  assuming  an  ideal  gas  mixture  and  the  equil. oration  of  translation¬ 
al  and  rotational  temperatures.  The  internal  energy  modes,  rotation,  vibration 
and  electronic  excitation,  are  considered  uncoupled;  and  a  rigid  rotator,  cut 
off  simple  harmonic  oscillator,  independent  cf  the  electronic  state,  is  assumed 
Excitation  of  vibrational  and  electronic  energies  are  treated  similarly  with 
terms  which  account  for  relaxation  and  chemical  reactions.  The  effects  of 
nonequilibrium  vibrational  and  electronic  states  on  chemical  rates  are  included 
in  the  coupling  analysis.  The  vibrational  relaxation  time  constants  were 
obtained  from  the,  Millikan  and  White  data  while  the  electronic  relaxation  time 
constants  were  determined  for  nitrogen  from  an  analysis  of  existing  shock  tube 
radiation  measurements.  The  computer  program  was  used  to  solve  for  the  non¬ 
equilibrium  flow  in  a  hypersonic  nozzle  and  for  eight  streamlines  in  the  invis- 
cid  flow  field  over  a  spherically  blunted  nine-degree  semiapex  angle  cone  at 
zero  angle  of  attack. 

Distribution  of  this  abstract  is  unlimited. 
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Ths  computer  program  described  in  thia  report  has  been  Witten  to  use 
only  the  cgs  system  of  units.  The  equations  for  vibrational  and  electronic 
relaxation  time  constants,  internal  to  the  program,  and  the  equations  used 
for  the  input  of  chemical  rate  constants  are  only  correct  as  written  in  the 
cga  system.  To  allow  the  program  user  to  work  with  either  the  U.3.  Custom* 
ary  System  (tf.S.)  or  the.  International  System  of  Units  (S.I.),  a  stbroutina— - 
wculd  have  to  be  written  to  convert  all  input  data  is^-cg5Vraans,,'  ':A  related 
subroutine  could  be  used  to  convert  cge. output  data  into  the  desired  system. 
The  units  previously  listed  would  furnish  the  basis  of  these  unit  conversions. 

Difficulty  may  be  experienced  in  converting  chemical  rate  constants,  ex- 
passed  as  a  .  unction  of  temperature  only,  from  Rankine  to  Ke7*vin.  Assume  - 
that  the  rate  constant  is  known  in  the  proper  form  based  bn  temperature  in  °R, 

fa  *-tc ooC 

K  *  AT*  «.  '"’T 


To  find  ths  input  -data  constants,  a,  b,  and  c,  -which Hill  be  used  within  the 
program  with  temperature  in  °K, 
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it  can  be  shown  that 

a  =  (1.8)BA 


c  =  C/1. 8 

a  further  correction  to  a  is  needed  depending  on.  the  units  used  for  species 
concentration.  ^ 
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SECTION  % 


PRODUCTION 


The  object! vp  of  this  study  was  to  modify  an  existing  ohd^dMehslonal 
reacting  gafj;  ;f ion  computer  program,  Reference  1  to  include  vibrational  and 
electronic  relaxation  processes  and  chemical  coupling.  This  older  program 
had  the-  capability  oi  predicting  test  section  conditions  following  the-npp- 
equilibriuni  exparoion  of  air  in  the  nozzle  of  an  arc-heated  experimental  hyper¬ 
sonic  facility;  It  could;  also  be  applied  to  the  shock layer  analysis  o£  boh-* ' 
equilibrium;  flow  over  test  models  placed  In  such' a  nohequilibrium  and  nbnuhi- 
form  test  section  flew  as  well  as  full  scale  atmospheric  entry  flow  fields.  '■ 
Nonequilibrium;  vibrational  and  electronic  excitation  can  bp' ejected:  to  become : 
more  importaiit  as  the  enthalpy  levels  of  test  facilities  axe  increased  id  ah 
effor1  to  moro  closely  simulate  the  full  scale  atmospheric  entrjr  environment  ? 
of  current  and  future  space  vehicles.  This  consideration  led  to  the  require-- ; 
merit  for  an  improved  reacting  gas  program.  J  '  '  - 

A  review  Of  the  current  state-of-the-art  showed  that  the  early  work  Of 
Tree nor  and  &rrone  o/ra  nonpreferential  coupled  vibration  model.  Reference '2  ,. 
had  been  continued  to  develop  a  simple  preferential  model.  Reference  3  .  In 
1965  'Premier  presented  a  sophisticated  analysis,  Reference  U  >  of  preferential 
dissociat  :cu  from  excited  vibrational  levels  which  are  not  in  a  Boltzmann  dis¬ 
tribution.  The  model  proposed1  in  Reference  3  was  adopted  for  this  -study  s‘ihce: 
it  was  OfjiiSidered  to  be,  sufficiently  accurate  for  engineering  analysis,  con¬ 
sidering  the  many  uncertainties  in  the  basic  data  needed  for  these  calculations, 
and  was  net  much  more  difficult  to  program  than  the  original  nonpreferential  "  ; 
model  of  Reference  2  No  foreseeable  difficulty  was  anticipated  in  extending 
the  analysis  to  include  electronic  processes.  V 

the  basic  design'  of  the  nonequilibrium,  program  Was-  to  remain  unchanged 
•eyth  though  it  completely  new  computer  program  was  formulated  and  written..  Sub¬ 
routines  would  be  added,  where  necessary,  to  -cover  the  new  logic  required  by 
the  inclusion  of  finite  vibrational  and  electronic  relaxation  processes.  /The 
basic  reasons  fur  writing  a  new  program  stem  from  the  need  to  expend  the  chew* 
leal,  model  to.  include  more  species  and  reactions  and  the  need  to  include  a 
greater  : .umber,  of  integrated  variables  Iruhhe  form  of  vibration  ai;.d  electronic 
parameters.,  Both  effects  Resulted  in  changes  to  the  size  of  storage  arrays 
and  the-  indexing  system  so  extensively  used  throughout  the  program. 

"  .'*<?  .  *  ^  i  >  *  -  t  * 


SECTION  II 

THEORETICAL  -FORMULATION 


1.  BASIC  APPROACH 

An  effort  was  made  in  the  new  program  to  formulate  this  input  of  cheDdca^",, 
species  in  a  manner  that  will  allow  the  chemical  model  to  be  easier  aitero4  ,t 
or  replaced.  In  this  way,  it  will  be  possible  to  separately  study  various  V 
chemical  models  such  as  an  air  (oxygen-nitrogen )  model,  a  hydrogeri-oxyger  ; 
model,  a  carbon-hydrogen-nitrogen-oxygen  model,  etc.  ,  it  wili  also  #o^de 
flexibility  in  performing  future  studies  of  the  relative  importance  of  various 
species  and  reactions  included  in  the  chemical  model.  At  present,  our  primary 
interest  lies  in  the  air  model  which  will  be  enlarged  from  its  present  state. 
Provision  is  made  within  the  program,  to  include  up  to  fifty  reactions  and" 
twenty  species.  The  two  new  species  which  are  added  to  the  air  chemistry,, 
model  are  0”  and  ($•  /  ; 

Along  with  the  enlargement,  of  the  chemical  model,  it  is  also  necessary 
to  make  provision  within  the  program  for  finite  changes  in  the  vibrational, 
energy  of  each  diatonic  species  and  changes  in  the  electronic  energy  for, all 
species,  both  monatomic  and  diatomic.  This  implies  a  step-wise  integration 
of  vibration  and  electronic  energies  similar  to  that  currently  being  performed 
on  the  chemical  species  concentrations .  The  inclusion  of  these  other  energy  y 
modes  in  the  formulation  further  complicates  the  flow  equation  which  relates 
the  changes  in  chemical  composition,  and  now  changes  in  Vibrational  and  elec¬ 
tronic  energies,  to  the  fluid  properties  through  a  differential  expression 
for  the  local  translational  temperature,  Tt  .  In  the  revised  version  of  the 
program,  vibrational  and  electronic  ..energies  are  associated  with  their  own 
individual  temperatures,  “firt  ,  and  jei  ,  respectively,  which  are  different 
from  the  translational  temperature.  The  relationships  between  these  new. . 
temperatures  and  the  translational  temperature  will  become  more  apparent  in  a 
later  development  of  the  flow  equation.  The  program  is  formulated  in  a  very 
general  way  so  that  various  assumptions  can  easily  be  made  as  to  the.  vibra¬ 
tional  and  electronic  state  of  the  model. 

Two  forms  of  coupling  vibrational  relaxation  to  chemical  relaxation  were 
included  in  Treanor's  model.  One  form  pertains  to  the  effect  of  vibration  on 
dissociation;  whereas  the  other  pertains  to  the  effect  of  dissociation  on 
vibration.  The  first  effect  is  introduced  in  the  form  of  a  coupling  factor, 
Viq,  (one  for  each  applicable  diatomic  species),  which  is  used  to  modify  the 
equilibrium  wue  of  the  cheiical  reaction  i^te  constant  for  those  reactions 
in  which  chemical  bonds  are  broken.  F*  -jsi  this  standpoint,  the  coupling,  factor 
may  be  applied  to  either  the  forward  or  reverse  reaction  rate  constant.  Al¬ 
though  the  effect  of  electronic  excitation  on  the  chemistry  processes  is  not 
ae  well  understood,  provision  is  also  Bade  in  the  program  for  including  the 
effect  of  electronic  relaxation  on  the  chemical  rate  constants  in  the  form  of 
a  similar  coupling  factor.  At  present,  all  of  these  electronic  coupling  fac¬ 
tors  will  be  set  to  unity,  shunting  out  their  effect,  because  of  a  present 
day  lack  of  knowledge.  However,  it  will  be  possible  with  this  type  of  formu¬ 
lation  to  maintain  complete  generality  within  the  program. 
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the  second  form  of  coupling  employed  by  Treanor  occurs,  through  ter ms  in¬ 
cluded  in  the  vibrational  rate  expression  foi*  each  diatomic  species*  these 
terms  pertain  to  the  vibrational  energy  lost,  ,  or* gained,  S ir^  ,  by  the 
destruction  or  formation  of  a  molecular  specie^, respectively.  Similar  terms 
expressing  the  gain  or  loss  of  electronic  energy  from  the  formation  and  des¬ 
truction  of  chemical  species  are  included  in  the  hew  formulation*,  the  method 
of  evaluating  such  terms  for  the  electronic  case  is  still  quite  questionable 
and  it  is  necessary  to  again  shunt  all  of  these  terms  from  the  program  until 
a  better  understanding  of  them  has  been  obtained.  , 

Another  part  of  the  program  modification  deals  with  the  expressions  for 
the  individual  species  thermodynamic  properties.,  enthalpy  and  specific  heat. 

In  the  old  version  of  the  nonequilibrium  chemistry  program  species  thermodynamic 
properties  were  based  on  an  equilibration  of  the  vibrational  and  electronic 
temperatures  to  the  local  translational  temperature.  The  new  expressions  are 
not  as  restrictive  and  contain  provision  for  including  individual  species  vibra¬ 
tional  and  electronic  thermal  effects.  The  new  expressions  are  based  op  a 
rigid  rotator  cutoff  harmonic  oscillator  model  and  it  is  presumed  that  there 
is  ho  coupling  between  the  various  internal  energy  modes. 

The  final  modification  is  the  inclusion  of  Treanor »s  integration  proced¬ 
ure  as  described  in  deference  5  .  This  procedure  is  designed  to  help  over¬ 

come  the  difficulty  associated  with  the  numerical  integration  of  a  set  of 
coupled  differential  equations  with  greatly  differing  time  constants.  In  this 
type  of  problem,  commonly  referred  to  as  the  “.stiff"  equation  problem,,  ,step  , 
size  is  determined  by  the  fastest  rate  and  the  region,  of  integration  is  set  by 
the  slowest  rate.  The  term  ''stiff'  refers  to  an  extreme  sensitivity  of  the 
rate  of  a  .parameter  to  small  deviations  of  the  integrated  value  from  the  true 
curve.  Treanor 's  procedure  is  a  modified  version  of  the  standard  fourth-order 
Runge-Kutta  integration  technique  which  is  specifically  designed*  to  minimize 
the  “stiffness"  effect  and  permit  a  larger  integration  step  oize.  In  regions 
of  the  integration  process  where  "stiffness "  is  not  a  problem,  Treanor '  s  phq^*, 
cedure  becomes  identical  to  the  Runge— Kutta  procedure,  Thus,  the  new  integra¬ 
tion  scheme  is.  designed  to  furnish  the  same  accuracy  but.  at an  Improved  program 

run  time. 
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This  section  presents  the  reformulation  of  the  theory  used  in  the  nonequil— 
ibrium  chemistry  program  to  include  finite  vibrational  and  electronic  relaxa¬ 
tion  processes'  and  the  Coupling  effects*  between  these  energy  modes  and  the 
chemistry.  We  begin  by  presenting  the  three  fluid  conservation  equations. 

Conservation  of  Mass 

pa, A  *  M  (l) 
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Conservation,  of.  Momantum 
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Conservation  of  JEhergy 
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To  these  equations  we  add  expressions  relating  to  the  thermodynamic  state  of 
the  "gas  media  in  the  form  of  the  state  equations  and  the  definition  for  on- \ 
thalpy  of  the  gaseous  media,  (y-. 
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Equation  of  State 
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Definition  of  Enthalpy 
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These  five  equations  are  taken  directly  from  Ref  erencel  with  the  exception 
of  the  modification  to  Equation  (5)  to  include  the  separate  effects,  of  the 


The  expression  for  the  chemical  rate  .of  change  of  a-speciesaole/maos 
ratio  with  distance  was  basically  defined  in  Reference  1  but  due-  to  the  lii-  ' 
elusion  .of  vibrational  and  "electronic  modes  of  excitement,  it  has  been  modified 
accordingly:  ‘  ~  ,  ...  .  /' 
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Chemical  Rate  Expression 
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The  coupling  effects  of  the  vibrational  ard  electronic  modes  on  the  chemistry 
appear  in  the  forward  and!  backward  chemical  reaction  rate  constants,  Kf^and 
These  coupling  effects  are  introduced  in  the  form  of  coupling  factors 
Vtr^  and  V/e*  for  the  vibrational  and  electronic  modes,  respectively.  The  net? 
rate  constants  are  defined  as  follows: 
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The  second  form  Of  coupling  discussed  by  Treanqr  in  References  2  and  3 
pertains  to  the  chemical  effect  on  vibrational  energy  rates  due  to  the  format 
tion  and  destruction  of  species.  The  relationship  for  the  vibrational  energy 
rate,  shown  in  Equation  (8)  below,  is  derived  from  an  energy  balance  as  des¬ 
cribed  in  Appendix  I.  For  a  given  mass  .of  gaseous  mixture,  .there  safe  three 
contributing  factors  to  the  rate  of  change,  of  the  vibrational  energy  state  of 
the  gas.  One  effect  is  due  entirely  to  molecular  excitement,  no  consideration 
being  given  to  chemical  reaction.  The  other  two  effects  account  for  the  gain 
and  loss  of  vibrational  energy  due  to  the  formation  and  destruction  of  chemical 
species,  respectively.  These  three  terms  appear  in  the  braced  quantity  in 
Equation  (8)  below.  The  remaining  terra  in  the  equation  pertaining  to  the.net 
chemical  change  is  derived  from  the  fact  that  the.  energy  rate  given  by  -Equation 
(8)  is  based  on  a  mass  of  species  "i*'  rather  than  a  mass  of , mixture-  basis-*  It 
accounts  for  the  fact  that  the  mass  of  species  'i"  changes  with  distance  along 
the  streamtube,  whereas  the  mass  of  mixture  remains  constant..  The  resultant 
vibrational  energy  change  with  streamtube  distance  is  expressed  by: 


Since  for  any  species  a  vibrational  temperature  is  defined  by  the  vibrational 
energy,  the  vibrational  temperature  derivative  is  computed  and  integrated.  This 
avoids  an  iterative  calculation  of  vibrational  temperature  from  the  vibrational 
energy.  Thus,  *  , 
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where  the  first  iera  on  tfa,  right  hand  side  is  computed  in  closed  fora  frcm 
Equation  (29)  .  c 

The  terms  £ sv  both  opcciesf ornaticn  and  destruction  are  posi¬ 

tive  quantities  in  Equation  (8)  and  are  defined  as  follows: 


Combining  Equations  (  9  )  and  (10),  as  shown  below*  yields  the  net  change  in 
chemical  species  with  streamtube  distance  which  is  identical  with  the  expres¬ 
sion  in  Equation  (6). 


Hie  relationship  shown  below,  and  similar  to  ihai  presented  in  Equation 
(8)  for  vibrational  energy,  expresses,  the  rate  of  change  of  electronic  energy 
with  streamtube  distance  which  is  used  in  the  computer  program.  The  form  of 
the  tern  dependent  on  the  electronic  relaxation  time  constant,  'te-  *  is  more 
fully  discussed  in  Section  II.6. 
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Since  for  any  species  an  electronic  temperature  i3  defined  by  the  electronic 
energy  the  electronic  temperature  derivative  is  confuted  and.. integrated  to 
avoid  the  iterative  calculation  tejnper&idre  from.*a  khovm.  ene^QE*. ,  ‘Shtts, 
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where  the  first  derivative  in  the  right  hand  side  is  computed,  in  closed.  form 
from  Equation  (14),  below.  Because  of  the  lack;  of  knowledge*  of  tgbj  for  other 
than  the  pure  species  Mfc  the,  computer  program  integrates  tie'  ^iec^ronic  temp¬ 
erature  of  N2  and  equates  all  other  species  electrode' temperatures  to  this 
value.  As  stated  earlier,  it  J.s  also  presently  beyond  the  state-of-the-art 
to  evaluate  the  terms  Se^apd  Therefore,,  they  will  be  set  equal  to 

j-h.  Equation  (12),  shunting  out  their  effect  and  leaving  for  the  .species 
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pie  energy  terms  Ee^Cl^and  E-s»-  are  obtained  from  the  folloiiing  expression  which 
is  based,  on  the  electronic  partition  function  evaluated  at  the  different;  energy 
levels  £  or  T*TVt  and  >  respectively:  -  .  *  '■  V/ 
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To  obtain  the  flow  equation,  which  relates  the  various  finitely  changing 
energy  modes  to  the  flow  .properties,  through,  the  energy  conservation  equation, 
one  begins  by  differentiating  the  enthalpjy  expression  presented  in  Equation 
(5).  This  yields  ' 


•  •>  ‘  ••  ft’ r ;  7s<  ■.  f. 

St  QdKMEi) 

J.%  £».  ds  £=,  aS 


Looking  for  the  moment  at  dkc Clt[ \Xvit ~Ttc  )  ,  we  find 


*•’  y  >•« 


i 

f 


ts a 


>  .I'. 

■4S  #S  dir  ■ 


Vhlch,  nhenanboUtuM  into  «*■**  ftjfc  ^  to  ^ 


oppression: 


f|  #  ^  i  *Ew«{  Jkvi 

«•"*.  AS  i*.|  v  -g~t 


+j^vik;a 


:s  '&>  (cm  4tit  s 

V*';,.  as 


*  X  >  * 


<fe.,/4s 

iSTlfg&Jta K"£  £*»«*<«  «>•  ®»  •*«*».  1.'  to 
*»«►  MSiUbrira  'Obtained  ..  - 

nay  be  streantuba  area,  aMiaatubi  .^OT^uST^a^ief1^  *?W  Wawter 
velocity  matching  relationship  ie  tto  eSt^  :**'  .  . 
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Tile  flow  equation  for  p 
through  the  use  of  the 
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Flow  Equation  ~  (Pressure  Hatch) 
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The  derivation  of  the  expression  for  area  match  is  more  involved  and;  reqUiteb 


eJA/ds. .  This  expression  is  then  substituted  ih  &iuktion  (l8k). , to;  produce  ,  '  , 
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One  more  differential  equation  relating  a  flow  property  to  ^he  anprn- 


reacting.  gaa  problem.  For  the  cases  of  velocity  an&  prqssyre  matching,  one 
can  use?  the  momentum  equation  and,  integrate  pressure  pr  yelocity,  respectively, 
as  shovm  below..  .  .  ’  •  -k! 
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In  the  area  match  case,  velocity  will  be  the  second  integrated  flew  parameter 
as  given  below:  ‘  V"1’ 
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This  section  deals  with  the  application  of  this  modified  theoretical  de¬ 
velopment  to  the  nbhequilibrium  problem.  The  first  task  is  to  define  a  set 
of  starting-  conditions.  In  the  case  of  external  flow  the  starting  point  is 
generally  assumed  to  lie  immediately  downstream  of  the  shock  front,  and  the 
pressure-match  relationship  for  ^Tv/ds  (Equation  (18b))  isused,  For  nozzle 
flow  Where  the  area-matqh  relationship  for  dTc/ds  (Equation.  (18c)  is  used,  , 
one  finds  that  a  removable  ^angularity  exists,  at  a  point  in  the  flow  where 
the  first  bracketed  quantity  in  Equation  (l8c)  goes  to  zero.  It  can  be  shown 
that  tills  bracketed- quantity  is  equal  to  ihe  difference  . between  the  frozen 
Mach  number  squared  and  unity.  Furthermore,  this  point  occurs  downstream-  of 
the  throat  and  it  is  necessary  to  start  the  nozzle  case  downstream  of  this 
singularity.  Figures  1  and  2  can  be  used,  to  determine  the  approximate  loca¬ 
tion  of  the  singular  point  in  the  nozzle  from  known  equilibrium  conditions.i 
There  is  no  such  singularity  in  the  pressure-match  relationship.  In  the  ex¬ 
ternal  flow  problem,  translational  temperature,  ,  is  based  on  immediate 
equilibration  of  the  translational  and  rotational  modes  of  excitement  behind 
the  -shock;  The  vibrational  temperatures ,  "lifi.  ,  for  all  diatomic  species  and 
the  electronic  tenQ>eratuPes,  Te^  ,  for- all  applicable  -species  'Wiil  be  set 
equal  to  the  free  stream  temperature,  fin  passing  through  the  shock  front  the 
chemistry  is  assumed  frozen  and  equal  to  the  free  stream  composition.  In  the 
nozzle  case  the  conposition  and  all  temperatures  will  be  set  equal  to  equilib¬ 
rium  values  for  tunnel  test  conditions,  '  \  , 


Having  described  the  starting  point  conditions,,  we  proceed  with  a  des~ 
cription  of  the  procedure  for  advancing  downstream.  This  procedure  is  des¬ 
cribed  below: 

.  .  •  •  ■  5  ■’* 


Step  1 


Step  2 


Knowing  the  initial  values  6f~ft  and  Tv^  ,  compute  the  vibra¬ 
tion-dissociation  coupling  factors,  V«r £  ,  for  each  diatomic" 


-species  from  Equations  (25),  (26),  and  (27). 
will  be  temporarily  set  equal  to  unity. 


The  terms  Ve£. 


The  forward  and  backward  chemical  reaction  rates  based  on 
translational  temperature  will  be  evaluated  for  each  reaction 
and  entered  into  Equations  (7a  and  7b)  along  with  the  terms 
computed  in  Step  1  above.  The  modified  rate  constants,  as 
functions  of  the  three  temperatures,  "Tt  9  lw~c  ,  andTci.  ,  will 
then  be  determined. 


Step . 3  -  Once  the  chemical  rate  constants  are  available,  one  can  proceed 
with  the  evaluation  of  the-  net  chemical,  rates  for  each  species, 
(<fcu/d&  )net,  from  Equation  (6)  and  the  species  formation  and 
destruction  terms  (4*4 /ds  Xit^a^  and  (dm/ As  ).h.Be^r  * 
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Figure  1.,  Equilibrium  Sound-  Speed  ifop  Mixture-  pif' Reacting 
62  and  0  arid  Inert  N2  (Reference  J.  Der,  NASA 
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Figure  2.  Ratio  of  Froaen  to  Equilibrium  Sound  Speeds 
For  Mixture  of  Reacting  O2  arid  6  and  Inert  % 
(Reference  J.  Der,  NASA  TR  R- 164) 


Step  3  —  (continued) 

from  Equations  (9)  said  (lo)  which  are  used  in  Equations  (8) 
and  (12)  for  the  vibrational  and  electronic  rates,  respectively. 

Step  U  -  The  vibrational  energy  terns  Errit0:  and&tr^  are  evaluated  by 
means  of  Equation  (29)  using  the  temperatures  It  and./Tv\;  , 
resisb'ct'ivelyv  /The  energy  t eras  ^5i  t  and^S^t  3  ar$  evaluated 
from  Equation  (29)  and  Tables^  VI  dm  VII.  The  vibrational 
relaxation  tics  ,Cfc*i)*i«  »  is  detemlned  from  Equations  (30)  and 
01).  Using  thlsand  -i^ormatloh  from  Step  3  above,  one  can 

_  evaluate  the jWbrationai  energy  rates  fticm  Equation  ($5  for 

,  each  diatimic  specie. 1  -  -  -  t.  ~  .  : 

Step  5  -  For  the  present  ,Se.^ ,  t  «hd  5^t:,  will  all  be  set  equal,  to 

each  other  and  the*  electronic  energy  rate  for  molecular  nitro¬ 
gen  will  be  computed  from  Equation;  (13)^  The  electronic  temp¬ 
erature  ofnitrogeh  will  be  integrated  and.  the  electronic 
temperature  of  all  other,  species  will  be  assumed  to  be  equal 
to  that  of  nitrogen.  The  terms  fife^CTr^ ) and  Se-  will  be 
evaluated  from  Equation  (14)  foi'Jl^^  ajjdTe,^  ,  respectively. 

Step  6  -  The  final  step  prior  to  evaluating  the  appropriate  flow  equation 

is  to  evaluate  *cpc  til)  and  ~tvi  ,  Tci)  for  each  species 

from  Equations  (22)  ami  (23).  % 

Step  7  -  Evaluate  the  appropriate  flow  equation  ,  Equation  Cls),  depend¬ 
ing  on  the  matching  parameter  used,.  5 

I  '  '  .  1 

Step  8  -  Equations  (6  )/ '  (^  );,  (12a),  (18)  and  one  of  the  Equations 

(19 —  21)  will  be  iorwardly  integrated  to  the  next  downstream 
point  in  the  streantube.  A  new  chemical  composition,  a  new 
set  of  vibrational  and  electronic  tejqoeratures,  a  new  trans¬ 
lational  temperature,  and  finally  a  hew  value  of  «.  or  p  will 
be  'detenainedatthe  new  station.  ;•  . 

»  '  *  5  1  * 
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Step  9  -  Since  the  matching  parameter  is  input  to  the  program,  its  hew 
value  can  readily  be  coeluted  at  the  new  downstream  point.  The 
remaining  flew  properties  at  the  new  point  are  obtained  from 
the  .mass,  conservation  equation  and.  the  equation  of  state. 

»  *,  %  .  ■  '  r  ■*»  j 

Having  performed  these  nine  steps,  .one  completely  defines  all  of  the 
required  parameters  at  the  new  station  .and  needs  rally  tc  repeat  the  procedure 
to  traverse  the  remainder  of  the  streanfcube.  ' 

3.  Tmaaminic  k®el 

‘Hie  expressions  for  the  thornodynaafe  properties  lvv  (“TJ  ,  Tr^ ,  Tc^  )  and 
for  the  uncoupled,  :utoff  harmonic  oscillator  model  are  ehram  below: 


In  the  expression  for  species  enthalpy  the  four  terms  within  the  braces  are  *  * 
translational  energy,  rotational  energy,  vibrational  energy,  and  electronic,  ^ 
energy*  The  term  h oL  represents  the  reference  or  heat  of  formation  energy  f°r 
the  particular  species.  The  vibrational  energy  is  based  ort  a,  cut-off  simple 
harmonic  oscillator  for  diatomic  sp8d.es.  to  ensure  consistency  with  theivibra~ 
tional  .coupling  analysis.  It  is  also  necessary,  3ince  noneqnilibrium; vibra-7 
tional  and  electronic  energies  are  being,  considered^ .-to  pepailate  the  yibrhtioil- 
al  and  electronic  energy  modes  from  the  translational  and  rotational  energy  . 
modes.  The  above  model  accomplishes,  this  separation.  Xt  assumes  a  rigid  ' 
rotator  with  rotational  temperature  eq’  to  the  translational  temperature.  v 
Diatomic  species  have  a  fiiced  number  o*  equally  spaced*  vibrational  levels,  l 
shown,  in  Table  I.  independent  of  either  translational  or  .electronic  tempera- \ 
ture.  These  data  were  obtained  from  Reference  '6  .  '<  ™  •  ~  ’ 

‘  ■  '  1'  •  i  ■» 
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,  ,  Table  I.  Vibrationhl  Energy  Constants' 1 


Specie 

:  M 

h 

°2 

NO  ,  ! 

N0+ 

W  i- 
■  0?”  .... 

.  *  '  . .  *  0  *  *  *% 

i  3374 

2256  . 

-  2719.  ■  .. 

■  339.7  : 

■  -  ,  3129  ■  .. ;. 

'  .  .  .2652  ■  ;  ! 

, .  1830  ,  '  ‘  > 

:  i  *  33  - 

'  "■  26* y 

27  : 

;  :  '  36  •  y 

:  - s,-  32 . V. 

28 

f  29  ‘  • ' 

The  electronic  energy  levels  are  assumed  as  cbnstants,  i^nd  ape  listed  in. 

Table  II,  These  are  based  on  date  presented  in  Ileferances  7  and  8  '  Energies 

hr,Vd  been  included  that  correspond  to  one  electron  excited  to  the  m,  -=  5  shell. 


The  higher  energy  levels  have  been  summed  into  grbups  using  energies  weighted 
by  their  degeneracies.  This  approach  follows  the  suggestion  of  Gilmore  who 
pointed  out  tha-  no  problems  with  series  convergence  will*be -encountered 
using  this  approach  while  good  accuracy  is  maintained  since  the  chemistry 
limits  the  number  of  molecule3  in  the  extreme  upper  levels*  A  comparison  ^d'th 
accepted  valued  showed  a  high  accuracy  of  species  enthalpy  computed  with  this 
model  at  550Q6K.  Subsequent  comparisons  have  been  made  for  an  air  mixture  up 
to  16,000°K.  An  examination  of  the  equilibrium  air  compositions  presented  in 
Reference  9  indicates  that  the  highest  concentrations  of  molecular  species 
at  high  temperatures  occur  at  the  highest  density  ratio  shown.  Thus, an  eval¬ 
uation  of  the  accuracy  of  computing  thermodynamic  properties  was  made  fop 
(?/po  ~  1»  where  f>0  is  air  density  at  1  atmosphere  of  pressure  at  6°C.  The 
percentage  error  ir  species  enthalpy  comp  -ed  with  the  above  model  was  found 
by  comparing  to  corresponding  data  in  Reference  7.  Mole  fractions  at  equilib¬ 
rium  for  the  species  N2,  O2,  and  NO  were  obtained  from  Reference  9,  The  error 
in  evaluating  the  enthalpy  of  an  air  mixture,  which  is  proportional  to  the 
product  of  the  species  mole  fraction  and  the  error  in  species  enthalpy  ,  was  com¬ 
puted  for  these  three  molecular  species.  The  results*  given  in  Table  II j, 
show  an  error  of  less  than  one  percent  which  decreases  with  increasing  temper¬ 
ature.  ; 


Table  III.  Evaluation  of  SD  Thermodynamic  Data 


k.  CHEMICAL  MODEL  FOR  AIR 

The  recent  publication  of  a  consensus  opinion  on  a  standard  chemical 
kinetics  model  for  flow  field  analysis.  Reference  10,  has  simplified  the 
problem  of  selecting  a  set  of  air  reactions  and  ihe  best  available  rate  data. 

It  has  been  decided  to  use  the  Bortner  model  for  air  as  presented  in  Tables 
IV  and  V.  The  computer  program  will  accept  up  to  20  species  and  50  reactions. 

The  sample  cases  have  been  run;  with  a  reduced  number  of  species  and  reactions, 
however  all  data  for  the  full  air  model  are  given  in  these  tables,  Bortner 
presents  two  sets  of  rate  data.  The  5 high  temperature'  data,  Table  XV,  are 
applicable  when  the  reaction  is  predominantly  endothermic,  while  the  'low 
temperature'  data,  Table  V,  are  applicable  when  the  reaction  is  exothermic. 

The  program  will  accept  either  set  of  data,  as  well  as  any  other  set  if  ex¬ 
pressed  in  the  same  form.  To  avoid  discontinuities  in  the  data  and  possible 
difficulties  in  the  integration  procedure,  no  provision  is  made  to  automatic¬ 
ally  switch  from  one  set  to  another  during  a  run.  Thus,  for  a  nozzle 
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expansion  flow  problem  where,  for  example,  recombination  reactions  become  pre¬ 
dominant  over  dissociation  reactions,,  th&  best  choice  is  the-  low  temperature 
set  of  data.  Fpr  an  external  flow  field  streamline  which  starts  at  a  shbck' 
wave  the  best  data  is  the  high  temperature  set,  -since  endothermic  reactions  • 
dominate  most  of  the  flow  field  and  Both  sets  are  equally  accurate-  near  -equil¬ 
ibrium.  Both  sets  of  data  are  consistent  with  the  same  equilibrium  cohstahtS  ,t 
so  that  identical  equilibrium  composition  will  result  from  use  of  either  hetV  " 
The  data  in  Reference  20  have  been  converted  into  values  which  are  consistent 
with  concentrations  expressed  in  g-moles/g^,.  Ths  rate  constants  are  given  by 
the  form  *  •  *  •  '  ’  -  ?r.r  ! 


cl  T  c. 


and  the  tables  list  values  for  a,  b,  and  c;.  -Due  to  the  fact  that  seme  values 
of  the  constant  a  are  exceedingly  large  and  may  present  machine  computational 
difficulties  in  the  evaluation  of  the  chemical  rate  constants,  all  values  of 
•'a'1  shown  in  the  table  have  been  reduced  by  a  factor  of  10°.  In  the  program 
the  rates  are  computed  using  the  exact  data  furnished  in  the  table  and  then 
multiplied  by  10°  to  obtain  the  correct  rate  constants.  -• 


5.  VIBRATIONAL  RELAXATION  EVALUATION  '  * '»  -  >  " 

In  this  subsection  the  form  of  the  expressions  for  the  coupling  factor 
and  the  energy  terms  will  be  presented  and  the  evaluation  of  the  expressions 
will  be  discussed  with  emphasis  oh  the  Bortner  chemical  model  for  air.  The 
vibrational  coupling  fadfcors yVvv*  have  been  thoroughly  defined  by  Treanor 
for  the  cut-off  harmonic  oscillator  model  in  terms  of  the  vibrational  parti¬ 
tion  function  evaluated  at  three  temperatures ,  "It ,  Tvr^  ,  and  an  effective 
temperature,  l?t  ,  given  by  the  relationship 
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Hie  vibrational  coupling  factor  is  obtained  from  the  expression 
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where  the  partition  function,  itself  is  defined  as. 


I  - 


K 


-  -  ---  v  r  —  *  ✓ 

-  -  -  -  -  *'  ,  _  /O  %  J1.-  ? 

"•  ■  ■  -  j  ^sC  v-'  f1' 

">  .■**'-  -  -  -  feA  . 

......  :y  .  „  !cf «  *v '  *>v¥i  ^  *• 

evaluated  at  the  appropriate  temperature.  The  negative  quantity  (-U)  can  he  '; 
considered  as  the  ‘’vibrational  temperature11  at  which  the  molecules  ar^-- formed 
by  recombination  in  Iteanor’s  preferential  dissociation  model.  .  For  the  limit-' 
ihg  case  of  U  .-«<©  ,  the  problem  reverts  to  the  nonpreferehtiai  model  ins 
Which  there  ip  an  equal  probability  <?f  dissociation,  frcja  any  vibrational  - 
energy  level.  ,  _  .  •..«•.<  .>*  ’*  “4  */•  - 

It  is  assumed  throughout  this  analysis  that  the  usual  eq\d.libriintt  rela*K  >' 
tionship  between  the  forward  and  backward  rate  constants,  as  defined  W|.bk,  ^  r 
still  applies.  Note,,  however,,  that  each  term  is  based  solely  on  the  trahsr  j 
lational  temperature.  -  .  .  .  ,  ) 


l<b^ 


*  ;  “‘vr>  -  ?.  '  “*  ?£V-/  / 


The  terms  ,  bv^  ,  and  ^0*  I  are  obtained  from  an  evaxuatich  pi 

the  energy  expression  for  the  cut-off  harmonic  oscillator  shown  below  u?idng 
the  temperatures  it  ,  ,  andl?^  ,  respectively..  .  ; -  /; 
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The  term^iQgt  in  Equation  (8)  is  also  obtained  from  Equation  (29)  by  ’  ,  . 
setting T»-U.  For  monatomic  species  ev-  -  <*> 

=  0,  and  VtT/  =1.  _ 

*  iw  V  ,  v. 

For  reactions  in  the  chemical  model  that  are  basically  not  of  thesis- 
sociation-fcecombinatipn  *ype,  it  may  be  necessary  to  evaluate  w 

by  ether  means,  than  Equation  (29).  The  preliminary  results  of  a  study  de^ 
signed  to  obtain  the  best  values  for  andGuv.:  are  presented  for  the 

Bortner  chemical  model  in  Tables  VI  and  VII  •»  Reactions  which  result  in  the 
net,  breaking  of  a  molecular  bond  are  presumed  to  involve  a  molecular  species 
having  an  and§u&,i  term  and  a  non-zero  or  .  The  deriva¬ 
tion  of  given  in  Reference  2  appears  valid  for  the  reaction  A2  +  X  + 

energy  jsS  2A  +  X  even  for  chemical  and  thermodynamic  nonequilibrium.  An  ex¬ 
ception  to  this  rule  is  made  for  reactions  in  which  a  molecular  bond  is  broken 
in  an  exothermic  reaction.  This  effect  is  observed  in,  for  example,  the 

reaction  . 

N  +  N  +  5.8ev^N^+  +  e~ 

Engineering  judgement  would  indicate  that  for  thes.  cases,  set 
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Reactions  in  which  there  is  ho  net  breaking  of  a  chemical  bond,  but  in  ■ 
which  an  atom  exchange  process  occurs,  are  stated  in  Reference  11  te  be  very 
effective  in  exciting  vibrational  energy.  The,  shuffle- reaptio^s,  for;  example,. 

,  ^  ^  .  .  •  »  -*£  H  ;  v  '0  VV  - 


are  believed  to  have  a  strong  effect  on  vibrational  excitation  of  the  species 
.•Njg  ,0£,  and  NO  .  To  account  lor  this,  effect jpn.  theseiaclecular  species  set 
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Reactions  in  which  a  -molecular  species;  gains  or  loses-  <an  electron  without 
a  net  breaking  of  a  molecular  bond  ’"are  presumed  to  have  no  vibrational  coup¬ 
ling  effects.  For  molecular  species  in  these  reactions.  Set 
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The  exponents  -fv&V  and  ,are  also  set,  equal  to  zero;  wider  any Jt£- 

the  following  additional  cases,  c.,  --  '  .  •' 

1.  The  specie  does  not  enter  into  the  reactipn  excspt  possibly  as  ^ 
a  catalytic  spscisa, 


,  jzy&ti 
.  ^  * 


t..  '  fee  .apociesie  mona tojnic'//. 

y.  It  is  desired' to  set  fee  coupling  factor  equal Jfco  1  without 
changing  ife  calculation  of  Vu^.  ,  - 
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The  present  a-st  of  exponents,  and 
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model  are  listed  below  for  all  susl^er  •»  values. 


for  the  Bortner  chemical 
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One  remaining  term  in  Equation  (8)  that  requires  defining  is  the  vibra¬ 
tional  relaxation  time,  Since  each,  vibrating^  species  can  conceivably 

collide  with  every  species  included  in  the  model,  it  is  necessary  to  find  a 
weighted  value  of  which  represents  the  entire  gaseous  mixture.  The  rela¬ 
tionship  for  this  weighted  value,  fr-  ,  -is  given  by: 
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The  derivation  of  this  expression  is  presented  in  Appendix  II  and  the  indi-  • 
vidual  relaxation  times ,T3rA , ftira  ,  etc.,  are  obtained  from  the  following 
equation,  taken  from  Reference  12. 
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where  is  the  total  pressure  for  the  entire  system  and  not  the  partial 
pressure  of  any  one  constituent, 

It  has  been  found  experimentally.  References  13 ,  and  14,  that  the  vibja- 
tional  relaxation  time  constant  for  expansion  flows  in  a  supersonic  nozzle  is 
much  less  than  is  predicted,  frcm  the  above  equation  which  correlates  measure - 
merits  of  the  relaxation  to  equilibrium  behind  shock  waves.  Since  the  explana¬ 
tion  of  this  effect  is  not  known  it  is  recosmnended  that  the  time  constants 
computed  for  the  above  equation  be  decreased  by  a  factor  of  15  for  nozzle 
expansion  type  flows. 

6.  ELECTRONIC  .RELAXATION  EVALUATION  .  \  '  '  .  .  '  ' 

The  electronic  energy  level  of  a  molecule  can  change  due  to  the  absorp¬ 
tion  or  emission  of  a  photon  or  by  various  Jypes  cf  inelastic  collisions  with 
other  particles:  neutral  molecules,  either  in  their  ground  states  or  in  ex¬ 
cited  states,  and  free  electrons.  The  rate  at  which  electronic  relaxation 
takes  place  depends  upon  the  numbers  of  the  different  species  present  and  the 
probabilities  of  each  type  of  inelastic  collision.  It  is  thus  coupled  to  all 
the  other  relaxation  processes,  in  principle. 


Attempts  to  describe  electronic  relaxation  have  lesorted,  of  necessity, 
to  various  assumptions  and  simplifications.  One  of  the  more  common  ones  is 
that  the  population  of  electronic  levels  is  Boltzmann-like. 

.  ■  s  ,  ,  <  ■> 
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The  electronic  temperature  "le.  can  be  regarded  simply  as  a  convenient  para- 
meter  with  which  to  describe  %he  distribution  of  electronic  energy. 

Allen,  et  al.  (Ref.  15}  measured  the  intensity  of  radiation,  I,  in  the. 
0.55  to  1.0  u,  region  during  the  passage  of  nitrogen  shocks  and  found  that 
I/I  eq  peaked  a  short  distance  behind  the  shock  and  then  decayed  to  unity 
(equilibrium).  The  parameter  (time  to  peak  intensity).  X  (initial  pressure 
of  unshocked  gas)  correlated  with  shock  velocity.  Time  to  peak  also  corre¬ 
lated.  with  the  vibrational  relaxation  time  for  nitrogen.  The  observed  radia¬ 
tion  behavior  was  compared  with  the  calculated  radiation,  from  the  nitrogen 
"first;  positive"  band  which  radiates  in  the  measured  wavelength  region. ‘  The 
calculations  assumed  that  %  was  equal  to  translational  temperature!  It  was 
found  that  observed  and  calculated  radiation  curves  agreed  during  the  later 
stages  of  relaxation,  but  riot  initially.  These  results  were  interpreted  as 
suggesting  that  electronic  temperature  remained  in  equilibrium  with  vibration¬ 
al  temperature  and  that  both  approached  the  translational  temperature  after 
the  radiation  peak.  . '■  . 

It  was  also  found  that  the  time  to  peak  radiation  obeyed  the  law  'tp.p  ~ 
*vh/. *  implies  that  the  electronic  relaxation  depends  ori  binary  col¬ 

lisions  (because  of  the  pressure  dependence)  and  that  the  effective  relaxation 
constant  is  a  function  of  the  translational  temperature,  at  least  in  the 

stages  (since  the  initial  translatiorial  temperature  is  a  function  of 
shock  speed).  This  kind  of  behavior  could  result  if  there  was  strong  coupling 
between  vibrational  and  electronic  energies,  since  the  vibratiorial  relaxation 
time  behaves  in  this  way  (Ref.  12).. 

-  </  k 

Later,  AHeri  (Ref.  16)  made  simultaneous  radiation  intensity  measurements 
on  shock  heated  nitrogen  in  different  wavelength  regions  with  a  resolution  of 
5A  and  determined  the  rotational,  vibrational,  arid  electronic  temperatures 
separately  by  using  nitrogen  itself  as  the .thermometric  substance.  The  rota¬ 
tional  structure  was  not  resolved  but  the  band  shape  for  the  N2+  "first  nega¬ 
tive"  system  was  used  to  deduce  7r  .  The  relative  band  intensities  of  th6 
same  system  were  used  to  obtain  Tv*  .  The  radiation  over  the  range  0.55  to 
l.lyx  was  used  to  deduce  the  electronic  temperature  from  the  N2  "first 
positive"  band  system.  The  results.,  for  a  shock  speed  of  6.4  mm/  U  sec,  show 
the  electronic  temperature  rising  more  rapidly  than  the-  vibrational  tempera¬ 
ture,  Te.  was  found  to  peak  (at  the  time  corresponding  to  peak  radiation 
ir  * ensity)  and  then  to  relax  downward  toward  the  equilibrium  temperature  with¬ 
out  ever  reaching  the  translational  temperature.  These  measurements  imply 
that  electronic  temperature  is  not  coupled  with  vibrational  temperature. 
Furthermore,  the. results  cannot  be  explained  in  terms  of  a  simple  relaxation 
equation  of  the  form,  ' 


-  J.. 


mUawa 
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Tzfjt)  =  relaxation  time,  assumed  a  function  of 
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This  is  because  the  electronic  temperature  peaks  and  begins  to  fall  while  ^ 
still  below  the  translational  temperature. 

It  should  be  eSgihasized  that  not  all  of  the  temperatures  were  determined 
from  the  same  species,.  ‘Piers  is  h,o  a  priori  reason  for  expecting  the  internal 
energy  modes  of  and  $2+  to  be  in  equilibrium'^  however  much  this  wbixld  ‘K 
simplify  the  situation.  '  '  •  ~  ^ ;  , P . 


Hansen  and  Chapin  (Ref,  17)  reach  similar  conclusions  from  ah' analysis' 
of  the  total  gas  cap  radiation  from  Small  spheres  fired  through  Sir 'in  a  bal¬ 
listic  range,  They  compared  experimental  values  of 

with  radiation  calculated  according  to  threat  different  assumptions:  '* 


•1. 


Te.  =  Tir 


Tile  population  of  excitedelectronic  states  rate-determined  by  t 
two-bedy  collisions,  with  neutrals  and  With  'electrons, ./the f^Utiwy  •“ 
collision  being  much  more' efficient.  f '  v  c*'\ ? 


For  the  third  type  of  calculation  it  was  necessary  to  introduce  rate'  constants 
for  electronic  excitation  which  wore  merely  plausible  estimates."  > 

It  was  found  that,  assumption  1  overestimates  ike  observed  radiation  by 
several  orders  of  magnitude,.  Assumption  2  gaye  closer  agreement ,  but  showed 
the  wrong  pressure  dependence.  Assumption  3  gave  reasonable  agreement  over 
the  range  M*--  lii,000  -  22,000  ft/sec,  p  =  1  -  100  non  Hg»,  The  calculated 
effective  electronic  temperatures,  (for  assumption  3)  were  differeht  for  dach 
species  but  gave  qualitatively  similar  histories,  namely,  electrphic  temper¬ 
atures  rising  more  rapidly  than  vibrational  temperature,  but  peaking  before 
reaching  the  , (higher )  trahslatipnai  temperature.  ’tV  ’  "  - 

Available,  wind  tunnel  information  on  electronic  relaxation  gives  a  'some¬ 
what  different  picture .  McGregor  and  3rewer  (Ref,  18)  carriedodt  simul¬ 
taneous  measurements  of  electronic  excitation  temperature  ’TwjJ  (as  determined 
spectroscopically  from  Equation  02))  and  free  electron  translational  temper¬ 
ature  Taf,  (as  measured  by  an  electrostatic  probe)  in  a  low  density  Argon  "  k 
free  jot.  Here,  both  le.  and.*Tli  were  out  of  equilibrium  throughout  the  entire 
flow  field,  it  Was  found  that  "Tc.  ahdT**,  agreed  within  35^  and  that  both  were 
a  factor  of  three  higher  than  the  bulk  temperature  of  the  neutral  Argon:  atoms. 


The  role  of  the  free  electron  iu  brought  but  by  the  worftbf  Rurle  and 
Russo  (Ref,  lit).  They  measured  the aodiir,  atoafixcitatibh.|ieJJi^®rature:  Tfjj^ 
in  Argon  and  1$  N2-Argon  during  expansion  in  a  wind  tunnel  nozzle.  They* 
found,  that  for  pure  Argon,  the,- 't^pratwyaTk.  fhozeik  at  a  value 

corresponding  to  the  (calculated  Argon  translational  temperature  at  the 'throat. 
When  1$  N2  was  addedl^  dropped  appreciably  and  agreed  With  the  calculated 
N2  vibrational  temperature,  TV*  . 

The  authors  of  Reference  14  explain  these  result#  as  follows;  The  N2 
vibrational  temperature  is  influenced  primarily  by  collisions  with  neutral 
molecules  and  relaxes  at  a  rate,  which  is  tirngs  faster  thsm  would  be ;.pre~ 
dieted  from  shock  tube  relaxation  data*  Vibrationally  excited'  N2  is  very 
effective  at  exchanging  energy  with  free  electrons  and  the  latter  take  up  & 
temperature  #}  Hr  ,  The  free  electrons,  in  turn,  are  very  effective  (in 
comparison  with  neutrals)  at  exchanging  energy  with  sbdium  atons  so  thgt^j, 
approaches  T«?  .  The  only  meaeured... quantity, „.of  _cour8e,..is  1*ja,  - 

The  significance  of  free  electrons  is  also  brought  out  by  Morgan  and 
Morrison  (Ref.  19).  They  make  estimates  <?f  the  relative  electronic  excitation 
efficiencies  of  free  .electron  and  neutral  collisions,  arid  fiM  that  the  former 
is  orders  cf  magnitude  higher..  Thpy„  conclude  that  imiediateiy‘:hehihd  th#.;;  '*  V 
shock  electronic:  excitation  and  Subsequent  ionization  take  place  by  neutral' 
collisions  only.  However.,  when  the  electron  fraction  rises'  to  a  value  of  the 
order  of  10*"^,  electron  collisions  become  the  dominant  mechanism  for  both  « 
excitation,  and  ionization.. 


The  explanations  given  in  the  references  cited  above  are  in  conflict  and2" 
do  not  provide- a  basis  for  the  rational  seiectich  of  an  electronic  relaxation 
model.  It  was  decided  to  approach  the  problem  by  a  reranajysis  of  the  nitrb~ 
gen  shock  data  from  Reference  15,  using  the  computer  program  developedunder  0 
this  contract  as,  a,  working  tool,  Nitrogen  was  chosen  because  of  its  rela¬ 
tively  simple  chemistry  and.  because  it  forms,  the  major  constituent  of; air . 

The  work  consisted  of  thrqe  ;parts.:-  f  ..  .  *  '  .  "  V  * 

1.  Computation  of  the  apparent  Te  variation  behind  nitrogen  shock 
waves,,  based  upon  a*  consistent  calculation  of  the  chemical  and  yibra-  c 
tional  rates.,  and  a  comparison  .between  experimental  and  theoretically 

v  calculated  .radiation,  histprieis.  ,  .  \ 

**  t  ,  *  -  t  , /*  •  .. 

2.  Formulation,  of  a  plausible  electronic  relaxation  model  consistent 
with  the  Te.  variation  found  in  part  1,  including  a  numerical  estimate 

, .  •  pf  therrate  .constant.  ‘ 

3.  Confirmation  of  the  model  by  recalculating  composition  and 
histories  using  the  relaxation  rate  expression  developed  in  part  2. 

The  shock  conditions  of  Reference  15  are  given  in  Table  VIII.  They  were 
calculated  from  measured  parameters  using  the  ideal  gas  approximation  with 
V  w  1.4  (nitrogen  behind  the  shock  is  assumed  to  have  Tr  -  It  ).  It 
was  further  assumed  for  simplicity  that  the  velocity  ratio  across  the  shock 
had  the  limiting  hypersonic  value. 1/6.  Initial  temperature  (not  stated  in 
Reference  15  but  presumably  room  temperature)  was  taken  to  be  300oK. 
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Table  Y1II.  Nitrogen  Shock  Conditions  (Before  Relaxation) 


Parameter 

Case  1 

~  Case  2 

Shock  speed  (tm/ jjl  sec) 

'  4.98 

5*34  t  ’ 

Initial  pressure  (mm) 

10,0 

3.00  * 

Temperature  behind  shock  (°K) 

12,000 

13,700 

pressure  behind  shock 
(dynes/cm^) 

3.19xl06 

1.09xL06 

Mass  flow  per  unit  area 
(gm/cnrsec) 

7.47 

2.4° 

Since  electronic  excitation  was  the  factor  to  be  found,  electronic  energy 
must  be  approximated  in  these  calculations.  This  was  done  by  assuming  «.  ' 
"Kr(  N2 ) .  Because  of  the  small  amount  of  energy  tied  up  as  electronic  energy 
for  the  conditions  in  question,  this  approximation  should  haye  little  influence 
on  the  composition  and  translational  temperature  history. 

A  noneqyilibrium  flow  field  calculation  was  carried  out  for  both  of  the 
shock  conditions  listed  in  Table  VHIusing  the  area  matching  option.  The 
calculation  yielded  translational  temperature  "fi:,  vibration  temperature  Tir, 
density  p  ,  and  the  concentrations  Of  Ng  ,  N2*  ,  N  ,  #  ,  and  e~  a3  a  func- 

tion  of  distance  behind  the  shock.  -  -  . 

**  *  '  ‘  -  .  k  ,  >) 

For  each  shock  speed,  theoretical  values  of  radiation  intensity  iii  the 
0.5  to  1.6  micron  region  were  calculated  at  a  series  of  locations  behind  the 
shock.  At  each  location,  several  values  of  intensity  V .re  Calculated,  based 
oh  parametric  variations  of  the  electronic  temperature  over  the  expected 
range.  The  radiation  calculations  were  carried  out  using  a  program  RADIAT, 
developed  by  North  American  Aviation  (Reference  20).  This  program  is  designed 
to  give  spectral  and  total  radiation  intensities  (integrated  over  a  finite 
wavelength  range)  for  gas  mixtures  as  a  function  of  temperature,  density,  and 
composition  of  the  radiating  species.  The  program  was  designed  to  compute 
radiation  for  equilibrium  gas  mixtures,  that  is,  where  Tfc  -  "fir  =  ~Tc-  f 
However,  the  above  temperatures  merely  determine  the  population  of  species  in 
each  electronic  and  vibrational  level.  Therefore,  emitted  radiation  intensity 
for  molecular  band  systems  will  be  correctly  evaluated  when  the  input  temper¬ 
ature  is  set  equal  to  Tr  =  Te  even  WhenTt  Te, 

Hansen  (Ref.  21)  has  shown  that  when  IV  the  spectral  intensity  from 
molecular  radiators  is  altered.  However,  ah  inspection  of  his  results  indi¬ 
cates  that  the  integrated  radiation  over  a  major  portion  of  any  vibrational- 
electronic  band  should  give  results  which  are  insensitive  to  15r. 

Figure  3  shows,  a  sample  output  from  RADIAT.  The  program  calculates 
separately  the  spectral  contributions  from  the  first  and  second  positive  bands 
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of  %  »  the  first  negative  band  of  %  }  the  'photoelectric  continuum  of 
nitrogen  atoms,  and  the  $ie.e-free, Radiation  from  eieetrons..  fn  all  cases  the 
first  positive  band  radiation  was  found  to  ptt  dominant  for  the  cases  analyzed. 

The  absolute  radiation  intensities  calculated  in  this  .  way  were  divided 
by  the  theoretical  radiation  at  equilibrium  (based  upon  extrapolated  compos¬ 
itions  and  translational  iemperaturai  to  obtsiri  a  series  of  theoretical  values 
of  I/I  eq  corresponding  to  various  distances  behind  the  shock  and  to  various 
assumed  values  of  T!e.  ’!•  '  -  X 

Experimental  values  of  the  relative  integrated  radiation  intensity 
I/I  eq  in  the  region  dominated  by  the  "first  positive"  band  system  are 
given1  in  Reference  15.  The  wavelength  ranges  are  slightly  different  in  the 
two  cases':  0.55  -  1.1  microns  for  =  4»98  ma/^  sec  and  0.55  -  1.0  microns 
for  -  5.34  atm/ju.  sec.  This  is  not,  expected  to  nuke  al significant  difference 
for  the  present  purpose.  ,•  Vf 

A  comparison  of  the  calculated  and  experimental  radiation  data  yielded 
a  series  of  Tt  values  for  which  the  two  agreed.  TPe  result’s,  for  the  lower 
shock  speed,,  are  given  in  Figure  4»  together  with  calculated  values  of  trans¬ 
lational  and  vibrational  temperatures.,  - 


The  results  were  somewhat  surprising.  The  calculated  vibrational  temper¬ 
ature/  rose  much  more  rapidly  than  the  electronic  temperature  and  peaked  much 
earlier  (not  shown  in  Figure  4).  IKe  electronic  temperature  continued  to 
rise  until  it  reached  the  translational  temperature  (or  perhaps  the  vibra¬ 
tional  temperature,  since  both  were  very  close  together  a£  that  point). 


Tentatively  the  following  form  of  the  relaxation  rate  law  was  assumed: 


a  4^5  =  (Ji)  ~  6re.  +  && 

^S.CTvv 


(34) 


where 


p) 


,  -s 

is  a  relaxation  time  characteristic  of 
neutral  collisions 


p  is  pressure  of  neutrals  (essentially'  total 

'  -measure) 


pill  ls  Measure,  of  electrons 


Pjj)  is  %  relaxation  time  characteristic  ol 
electron  collisions 


Such-  an  expr  sssion  would  describe  an  electronic  relaxation  which  tracks  the 
translational  temperature  during  the  early  stages  when  electron  concentration 
is,  very  low  and  then  tracks  vibrational  temperature  in  the  later  stages  be¬ 
cause  of  the  electron's  much  greater  ^.citation  cross  section.  Since  electron 
concentration  was  very  low  ('—3  ?c  10"  ^  mole  fraction)  it  the  low  speed  case,, 
it  was  decided  to  attempt  a  fit  using  only  the  first  totm  of  Equation  (34). 


Even  with  this  simplification,  there  were  further  difficulties.  The  empijd- 
cal  radiation  data  (and  therefore  the  empirical.  Tg,  curves)  vero  not  resolved, 
for  very  short  distances  behind  the  shock.  Also,  the  relaxation  time  con¬ 
stant  must  he  regarded  as  an  unknown  function  of  quantities  which  vary 
throughout  the  chock. 


These  difficulties  ware  met  by  assuming  first. that  the  electronic  time 
constant  was  inversely  proportional  to  pressure,  (as  is  the  vibrational  time 
constant)  and  second,,  that  the  product  p'tfe  was  a  rather  weak  function  of  Tt 
which  could  be  adequately  represented  by  an  average  value  over  a  range  of  a 
few  thousand  degrees.  With  these  approximations  Equation  (34).  could°be 
written  .  -  '  -  0 


o 


and  integrated  from  the  first. point  at  which  empirically  data  wa3  obtained 
up  to  peak:  *-  >  .  l  v  ' 


.Shr*. 


^  ’ ; iikpg$ /  &•  ”  6el * 


(36) 


where 


Max  is  the  electronic  energy  at  peak 

>  > 

itffa  is  the  initial  electronic  energy 


The  numerical  integration  implied  by  Equation  (36)  was  carried  out  for 
both  shock  cases  listed  in  Table  YIII  ,  yielding  average  values  at  affec¬ 
tive  translation  temperatures  of  7480°K  and  9970°K,  respectively.  They  are 
w  ~m  in  Figure  5,  together  with  the  corresponding  vibrational  ties  constant 
fov  N2  as  calculated  by  Equation  (31).  Although  there  was  little  basis  for 
an  adequate  estimate  of  temperature  dependence,  the  results  wars  fitted  to  an 
equation  Similar  to  the  vibration  equation  in  fora  to  yield  j 

*  '  * 

Equation  (3.7)  was  then  used  in  a  recalculation  of  the  lower  speed  shock 
case.  Again,  the  results  were  surprising.  It  was  discovered  that  Equation 
(37)  grossly  overestimated  the  initial  rate-of-rise  of  the  electronic  temper¬ 
ature.  The  behavior  near  peak  was  qualitatively  correct  but  peak  temperature 
was  displaced  upward  and  closer  to  the  shock  front. 


It  was  concluded  that  even  though  the  electron  concentrations  were  very 
1,w»  the  second  term  i .  Equation  (34)  was  important.  The  relaxation  rate 
data  were  rc-interpretec*  in  this  way  by  equating  rates  at  the  peak  electronic 
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V  v  - 

tl.  i&z  fee  CTv)  —  fee  ’  ,{38>. 

teCTir,  p.4 


where  ■“ 

* 

't'ep^  -  Z,Tk*jP  «<p <[.ioS«V^3  (39) 


The  results  of  a  final  trial  calculation  using  Equation  (3&)  are  also 
shown  in  Figure  4  ("Tv*  and  Fc  rsmainsd  essentially  unchanged  for  ail  trial 
calculations).  The  .results  agree  very  closely  with  empirical  electronic 
temperature.  This  result  appears  consistent  with  the  interpretations  of 
Reference  14  that  electronic  temperature,  is  excited  by  the  action  of  elec¬ 
trons  and  that  the  potential  causing  relaxation  is  governed  by  the  vibrational 
temperature  of  nitrogen  molecules.  - 

It  should  be  recognized  that  the  results  of  this  part  of  the  study  are 
very  limited.  They  apply  to  only  a ,  small  range  of  shock  speeds  and  for  one 
type  of  chemistry,  Moreover,  the  suggested  temperature  dependence  is  very 
speculative  at  this  time.  Much  more  work  is  required  to  extend  the  model 
and  the  associated  raue  data  to  the  ranges  of  temperature  and  composition 
needed  for  hypersonic  flight  and  wind  tunnel  test  interpretation. 
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The  oldvepsion  of  the-  nonequilibrira  program  employed  a  fourth  order 
Runge-Xutta  integration  scheme  Rodified  fpr  the  Gill  correction.  Treanor, 
Reference  5  ,  has  indicated  that  this'  scheme  can  be  uneconomical  to  Use  in 
regionh  of  the  flow  field  where  the  rate  of  a  parameter  in  a  higher  coupled 
ayatea  is  extremely  sensitive  to  slight  displacements  of  its  in«sgr>t8d  value, 
fro®  the  true  curve.  ^Reference  5  ,  Treanop  defines  a'msihcd  for  Vll«?iati^. 
this  sensitive  situation  -and  permitting  a  larger  step  size  $0  be  uae<L!  in  the. 
integration  process.  When  the  rates  are  not  sensitive  to  the  local  values, 
the  new  integration  schema  becomes  equivalent  to  the  original  Runge-Kutta 
procedure.-'  '  \  «  ;  ’*'*  '  \  <  " 

A  brief  description,  summarising  the  calculation  required  in  Trcanor’s 
method,  is  furnished  below.  The  process  of-  integrating  from  one  station  to 
another  in  the  flow  field  is  carried  out  in  four  basic  steps.  The  integration 
interval,  H  ,  is  divided  in  half  and  the  four  basiu  steps  pertain  to  two  calcu¬ 
lations  at  each  half  step  as  illustrated  in  Figure  6  below. 

o'  <  "  *“  -v  v ^  ^  '  3fc  -  -  '  -  -  '» 


Safe- 

-Xwt  Jr 


x  -  distance,  along  stream¬ 
line  (independent  variable)' 

y  =  flow  parameter,  specie 
concentration,  vibrational 
energy  (dependent  variable) 

❖  / 

f  -  derivative  of  y  with 
respect  to  x,  dy/dx 


Figure  6.'  -Illustration  of  Ruhge-Kutta 
“  :  Integration  ScheMe  v 


First  Cycle:.  (First  Half  Interval) 


Step  1 


-  <  Knowing  at, ,  ,  and  ■£,  ,  confute  SL^ahd  vj^from 


tt*  •= 


+  k 


K  f +  it 


Steg_2  -  Using  at-^and  ,  evaluate  The  method  of  evaluating  the 
derivatives  was  described  in  an  earlier  section. 
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Second  Cycle  (First  Half  Interval) 


Stop  1  -  Knowing  eft,  ,  ,.and  compute  3^  and  t^from 

4=5  ^  '  _ 

H,-  -  •*,  +  $**  *  -  •■-  - 
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Step  2  -  Using  and  ^  ,  evaluate  :% 
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Step  4  r  Determine  Fj  arid from  the  relationship:: 
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First  Cycle;  (Second  Half  Interval) 
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Step  1  -  Using  h,  4;,.  ,  fi  ,  ,  P  Pj  >.  and  ,  determine  jL  from 

H  .  n  '  s,  -  "  '  i*  .  5*W  -■ ■  ;  -  v  J 
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Step  2  v-  Determine  tft^from 

■  .  aL4  “  n  ,  •  '■.  ’ 

Step  3  .  -  Use  V^,  and  tj.^to  evaluate 


Second  Cycle 


Half interval) 


Step  1  ■-  Ccmputs  %  from  the  relationship 

*'  W'  '  '  ,  ' 

c=  ..  Pi.  4:  A  ■  , ' 


Step  2  ~  Compute.  from  the  relationship:  « 

*  (£*.  4r  *[($i  *■  p| *i )  r  ttf *-*-  p£i,) 

-  (f%  4~  P^l  +*  (  f<,*+  Pjjj 
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Step  2  ~  Compute  o^=  at  +  K  .  (  ,  „■ 

This  completes  the  integration  cycle  for  the  currant  internal  and  the  process 
is  repeated  for  each  succeeding  interval  downstream.  Note  that  point  5  is 
identical  to  point  1  for  the  next  interval.  ‘ 

The  integration  step  size  is  determined  in  part  from  the  "stiffness” 
parameter, P ,  as  defined  in  Step  3  for  the  second  cycle  in  the  first  half 
interval.  Treanor  states  that  the  integration  procedure  will  be  stable  for 
positive  values  of  P  so  long  as  the  value  of  Pn  •3l£  75.  For  negative  P  the 
Treanor  scheme  is  inferior  to  the  Runge-Kutta  method  and  the  program  should 
automatically  resort  to  the  Runge-Kutta  method  using  an  appropriately  reduced 
step  size  criteria  of  Ph  5.  it  has  been  found  that  under  certain  condi¬ 
tions  step  size  js  greatly  overestimated  when  it  is  based  solely  on  the  stiff¬ 
ness  parameter.  In  addition,  it  has  been  found  possible  for  the  intermediate 
values  at  points  2,  3,  and  4  to  be  negative.  To  ensure  a  foolproof  integration 
a  second  test  is  made  on  the  parameter 

f  ... 

Local  values)  of  ir  and  <4.  are  used  at  points  1,  2,  and  3,  This  test  is  used  to 
set  the  original  step  size  estimate  using  values,  at  .point  1?  Step 'size  is 
selected  to  set  the  value  of  the  parameter  exactly  to  the  nominal  value  as 
listed  in  Table  DC  provided  the  step  size  does  not  increase  more  than  a  limit¬ 
ing  value.  The  following  teats  are  allowed  to  only  reduce  step  size  as  irfliy 
cated.  No  increase  is  permitted.  A  second  test  of  the  parameter  is:  made  as 
soon  as  properties  and  derivatives  are  found  at  point  2,  where  only  an  upper 
bound  is  used.  Step  size  is  reduced  by  successive  corrections  based  on  16% 
of  the  indicated  step  size  error  (using  the  nominal  values)  until  the  test  is 
passed.  The  integration  now  computes  properties  and  derivatives  at  point  3. 

The  Ph  test  is  first  applied  at  this  point,  using  50%  of  the  indicated  step 
size  correction  fusing  the  nominal  value)  to  pl2.ce  Pti  below  the  upper  bound. 
Prior  to  making  the  Pin  test  a  check  is  made  of  all  parameters  to  see  that  ifyi. 
and  (the  terms  in  the  denominator  of  the  Ph  test)  dp  not  differ  by  less 
than  30  in  the  seventh  and  eighth  significant  figures.  If  the  difference  is 
less  than  30,  that  parameter  is  dropped  temporarily  from  the  Ph  test  because 
of  the  good  possibility  that  it  may  prove  highly  inaccurate  and  produce  an 
erroneous  reading  of  integration  instability.  The 

ib 

y 

parameter  is  applied  as  before  to  possibly  further  reduce  step  size.  If  the,' 
upper  bound  is  not  exceeded,  the  integration  for  this  interval  is  completed 
until  values  at  point  5  are  found.  These  values  become  the, initial  values 
for  the  next  integration  interval  and  the  procedure  is  repeated. 

It  has  been  found  desixeable  to  put.  a  limit  on  the  increase  between  step 
sizes  for  adjacent  intervals.  The  computer  program  now  uses  a  limiting  in¬ 
crease  of  10%,  however,  the  RAM  B3  sample  cases  used  a  limiting  increase  of 
100%.  The  RAM  B3  sample  cases  were  also  run  with  the  upper  bound  and  nominal 
values  of  the  Runge-Kutta  Ph  parameter  set  to  10.0  and  8,0.  Experience  with 
the  nozzle  flow  case  proved  the  desireability  of  using  the  lower  values  to 
reduce  the  variation  of  step  size  over  a  number  of  integrations. 
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SECTION  XU  . 

COMPUTER  mO®m  DEVELOPMENT- 


'  -fr  -  *  '  *\v' 

This  section  describes  the?  operating  bhahacteristi-cs  of  the  program; 
presents  a  flow  diagram  of  the  main  program;  and  describes  the  program  input 
and  output.  An  input  storage  location  ;map  is  ;aiso  included, to  aid  in  the  i 
preparation  of  input  data  sheets.  -  .  „  „  '  -  >  *  ; 


1.  OPERATING  ENVIRONMENT 


The  program  was  developed  and  checked  out  on  the  IBM  70%  computer  sys¬ 
tem  at  tht.  Space  Division  Computing  Center.  The  program  was  written  in 
FORTRAN  IV  language  and  compiled  with  IBM's  FORTRAN  IV  Compiler  Program.  A; 
trial  run  with  the  program  has  been  made  on  the  existing  IBM  7044/70%  direct- 
coupled  digital  computer  system  at  Wright-Patterson  AF3,  Ohio,  and  the  program: 
proved  compatible  with  their  computer  system  with  one  exception,.  All  unused 
spaces  in  the  data  field  or  input  cards  must  either  be  filled  in  with  zeros 
or  a  decimal  point  must  be  used  in  order  for  the  read  format  statement  to 
operate  properly,  Examples  of  this  can  be  seen  in  the  sample  data  sheets 
shown  later  in  the  text.  ■  ”  . 

2.  PROGRAM  DESCRIPTION 

The  computer  program  consists  of  a  main,  program  and  fourteen  subprograms. 
In  addition,  certain  library  routines  are  also  employed  by  the  program. 


A  rief  descriptive  flow  chart  of  the  main  program  is  presented  in  Figure 
7  to  illustrate  the  basic  calling  sequence  of  the  subprograms.  One  can  ob¬ 
serve  from  Figure  7  that  the  program  contains  an  initialization  procedure  and 
a  cyclic  process  which  consists  of  evaluating  derivatives  of  certain  flow 
field  parameters  and  integrating  these  derivatives  forwardly  along  a  predeter¬ 
mined  path  to  obtain  a  new  set  of  properties .  The  completion  of  an  integra¬ 
tion  interval  requires  cycling  through  the  loop  four  times .  On  the  fourth 
time  through  the  loop  the  value  of  MK  is  altered  from  one  to  two  and  the  pro¬ 
gram  prints  the  flow  field  properties  at  the  end  of  the  integration  interval, 
The  program  then  checks  to  see  If  it  should  continue  integrating  or  stop  alt 
its  present  location.  The  program  is  set  up  to  evaluate  the  derivatives  of 
all  parameters  for  the  start,  of  the  next  integration  interval  on  the  fourth, 
time  through  the  loop.  Thus,  if  it  is  to  continue  integrating,  it  heed  only- 
return  to  the  beginning  of  the  loop  and  repeat  the  process.  The  initializa¬ 
tion  procedure  at  the  start  of  the  program  consists  of  reading-in  all  input 
data,  setting  up  certain  working  arrays,  and  computing  derivatives  of  the 
integrated  parameters  at  the  starting  point.  The  selection  of  a  starting 
point  will  be  discussed  later  in  more  detail. 


'ltie  purpose  of  each  subprogram  is  briefly  discussed  below.  The  sub¬ 
program  BEGIN  is  responsible  for  reading-in  aii  input  data,  establishing 
certain  working  tables  and  arrays,  and  writing, -out  the  input  data  and 
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initial  starting  properties.  The  function  of  subprog: -am  THERMO  is  to  compute  , 
the  thermodynamic  properties,  enthalpy  and  specific  heat,,  for  each  chemical 
species  as  functions  of  the  three  temperatures  Tt,TIr.  4  andTe^  .  In  addition, 
it  evaluates  the  forward^  and  backward  chemical  rate  constants  based  solely  on 
translational  temperature.  Subprogram  VIB  is  basically  concerned  with  com¬ 
puting  individual  terms  U3ed  in.  evaluating  the  vibrational  energy  rates  for 
the  molecularjjpecies.  The  terms  computed  in  this  subprogram  are  e©  , 
£TiQ  >C*i  ,  ,Qit£  ,  and  ,  The  subprogram  ELECT  is  responsible  for 
obtaining  similar  terms  ( ,  and  te  )  used  In  the  electronic  rate 
expression  for  nitrogen.  Subprogram  CHMRAT  has  several  functions..  Firsts  it 
applies  the  vibrational  and  electronic  coupling  factors  to  the  forward  apd 
backward  chemical  rate  constants  computed  in  subprogram  THERMO;  These  rate 
constants  then  become  functions  ofTfr,TVt  ,  and .  .Second,  it  obtains  the 
chemical  formation  and  destruction  terms  for  each  species  and.  at  the  same  Jtiine 
combines  these  terms  with  the  ©r*_.,  ££vr*  ,  arid-^i/^  terms  used  in  the  vibra¬ 
tional  rate  expression,.  CHMRAT  also  combines,  the  chejaical  formation  and  ^des¬ 
truction  terms  used  to  form, species  net  chemical  pates  which  are  integrated  v 
parameters  in  the  program.  Subprogram  VIBRAT  evaluates  the  vibrdtional  energy 
rate,  dhir^  /is  ,  and  the  corresponding,  vibrational  temperature  rate,  /tts  , , 
for  each,  applicable  molecular  species.  Subprogram  ELCRAT  is  similarto  VIBRAT 
in  the  sense  that  it  evaluates  the  electronic  energy  rate,  d  ft*  /a  ,  and  the 
derivative  of  electronic  temperature,  die  /i$  .  The  program  is  temporarily  re¬ 
stricted  to  calculating  the  electronic  derivatives  for  only  the  species,  fe. 

The  last  subprogram  to  be  called  in  the  derivative  evaluation  sequence  is  the 
subprogram  PRPDER.  This  routine  uses  the  previously  calculated  i  vibrational, 
electronic,  and  chemical  rates  in  a  general  flow  relationship  for  dlb /ik, ,  One 
final  derivative,  d«.  /ds ,  is  also  evaluated  in  the  subprogram  PRPDER;.  Both  of 
these  derivatives  are  functions  of  the  snatching  parameter,  which  may  be  stream- 
tube  area  or  the  pressure  distribution  along  the  streamtube.  Subprogram 
ML00K  obtains  the,  value  of  the  matching  parameter  as  a  function  of  the  local 
position  along  the  streamtube.  Sinee  s  is  the  independent  variable  in  the 
analysis,  iocal  streamtube  coordinates  X  and  Y  are  obtained  from  the  newly 
integrated  value  of  s  in  this  subprogram.  The  subprogram  DERVP  is  called  by 
PRPDER  to  obtain  the  derivative  of  the  matching  parameter.  It  use3  a 
Lagrangian  three-point  fit  to  evaluate  the  derivative. ' 


The  final  subprogram  called  is  the  function  MARK.  This  subprogram  is 
responsible  for  integrating  the  previously  evaluated  derivatives  using  a 
fourth-order  Runge-Kutta  integration  procedure.  It  also  contains  an  elaborate 
testing  scheme,  which  was  discussed  in  detail  in  a  previous  section  of  this 
report,  for  altering  step-size  while  maintaining  accuracy  and  control  over 
the  stability  of  the  calculations.  There  are  four  branches  within  this  sub¬ 
program;  one  for  each  portion-  of  the  half-cycle  discussed  ift  the  theoretical 
section  on  the  integration  procedure.  On  the  fourth  time  through  the  sub¬ 
program  during  each  integration  step,  MK  is  changed  from  one  to  two  signaling 
the  completion  of  the  integration  procedure  for  that  step  ahd  allowing  the 
main  program  to  be  routed  through  the  print  -statement*  ’  , 

3*  OPERATING  INSTRUCTIONS  .* 

There  are  several  items  which  are  pertinent  to  the  proper  operation  of 
the  program.  For  instance,  the  value  of  XSTOP,  which  is  the  point  at  which 
the  program  stops  integrating,  must  lie  within  the  tabulated  matching 
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parameter  data.  It  is  also  necessary  to  make  the  input  valued,. of  the  coordi¬ 
nates  X  and  I  equal  to  first  values  of  TABX  and  TAB!  and  the  flow  property 
(pressure  or  area)  equal  to  the  first  value  of  the  matching  parameter, ='TABAP. 
In  this  way,  streamtube  distance,  s  ,  always  starts  . at  zero,  •  v 

’  '  '*  Ss*  .  '  4-'  * A 

If  the  print  indicator  LPRNT.  is'  set  to  zero,  a  limited  printout  of- input 
data  is  obtained.  This  basically  consists  of  input  data  that  may  change,  from 
case  to  case  when  the  chemical  model  remains  fixed.  For  example, -each  new 
case  will  have  new  starting  point  data  and  a  new  description  of  the  matching 
point  parameter  with  streamtube  distance.  On  the  other  hand  if  LPHNT  —  1,  a 
full  printout  of  all  input  data  is  obtained.  Two  other  indicators  that  re¬ 
quire  an  explanation  are  IND  and  UV„  If  IND  =  0,  streamtube  area  is  used  as 
the  matching  parameter.  For  a  pressure,  match,.  IND  is  set  equal  to  one.  rtThe 
indicator  UV  is  used  to  signal  the  program  to  follow  either  the, -prefer ential 
or  nonpreferential  dissociation  model  derived  by  Treanor.  The.  term  nonprefer¬ 
ential  applies  to  the  fact  that  dissociation  may  occur  .ap.  any  energy  level;..- 
whereas  the  term  preferential  applies  to  the  case  where  dissociation  occurs 
to  a  greater  extent  at  the  higher  energy  levels.  The  nonpreferential  model 
is  obtained  when  UV  =  6, ,  •  .  .  • 

The  program  is  presently  set  up  to  handle  a  limited  air  model  chemistry 
in.  line  with  the  current  AFFDL  tunnel  requirements.  The  species  present  in. 
the  model  are  N2,  0 g,  NO,  NO**  N>  Qf  e“,  and  Ar.  The  reactions;  included  in 
the  current  model,  with  reference  to  Tables.  IV  and  V,  are  1,.  Zy  3,  4„.  5,  6, 

7,  14,  15,  and  20.  Obviously,  one  must  always  insure  that,  there  is  compati¬ 
bility  between  the  species  and  reactions  used  in  the  program.  The  program 
has  the  provision  for  eliminating  both  bpecies  and  reactions  so,  long  as  this 
compatibility  restriction,  is  observed*  Species  are  eliminated  by  setting 
their  initial  concentrations  to  zero;  uid  reactions  are,  eliminated  by  setting 
their  forward  stoichiometric  coefficients  to  zero.  Bor.tner  gives  two  sets 
of  chemical  pate  data,  one  for  high  temperature  use  where  the  reactions  are » 
basically  endothermic  and  the.  other  fort  low  temperature  use  where  the  reac¬ 
tions  are  basically  exothermic,.  The  data  for  the  program,  parameters  FA,  EB, 
FC,  BA,  BB,  BC,  and  XMU,  all  pertaining  to  the  evaluation  ..of  the  cheMcal  - 
rates,  must  be  changed,  depending  on  which  case  is  involved. 

There  is  a  provision  within  the  program  for  adjusting  the  vibrational 
time  constant,  TAUV.,  when  a  nozzle  flow  case  is  being  considered,  •  As  stated 
ir.  the  theoretical  section  of  this  report,  . the  vibrational,  time  constant  for 
nozzle  flow  has  been  observed  to  be  about  fifteen  times  smaller  than  that  for 
external  flow  field  studies.  The  input  constant  TAUADJ  is  used  to  adjust  the 
time  constant  accordingly,  i  _ 

Initial  starting  conditions  for  the  program  will  either  be  representa¬ 
tive  of  a  point  immediately  behind  the  shock  in  the  case  of  external  stream¬ 
line  flow  or  a  point  slighvly  downstream  pi  the  throat  for  the  nozzle  case. 

For  the  shock  case,  one  assumes  that  translational  temperature  is  immediately 
equilibrated  and  that  the  vibration  and  electronic  temperatures  remain  equal 
to  the  free  stream  temperature.  It  is  also  assumed  that  chemical  dissociation 
does  not  occur  in.  passing  through  the  shock*  Therefore',  the  chemical  compos¬ 
ition  is  the  same  as  for  the  free  stream  air*.  The  remaining  flow  conditions 
are  then  established.  For  the  nozzle  case,  the  flow  is  assumed  to  bp  in 
thermal  and  chemical  equilibrium  and  Mfrozen  ^  T.‘ 


o 


4  .  INPUT  DATA  t;  -  - 

The  'Input  data  to  the  program  are  entered  according  to  the  storage-  loca¬ 
tion  nap  presented  in  Table  X.  All  parameters  identified  with  a  single  aster¬ 
isk  are  basic  input  information  to  the  program  and  generally  remain  constant 
except  when  modifications  are  made  to  the  chemical  model.  The  parameters 
marked  with  a  double  asterisk  are  .also;  input  and  basically  change  from  case 
to  case.  The  FORTRAN  symbols  appearing  in  Table  X  are  defined  in  the  list  of 
Abbreviations  and  Symbols  at  the  beginning  of  this  report,! 

A  complete  set  of  sample  data  sheets  for  the  RAM1B3  streamline  No.  5  is 
presented  in  Appendix  III.  Data  sheets  for  the  other  streamline  cases  are 
presented  in  Appendixes  IV-X,  as  changes  to  the  input  data  for  the  stream¬ 
line  No.  5  case.  The  input  data  sheets  for  the  nozzle  case  are  presented;  in 
Appendix  XI  and  differ  from  the  RAM-B3  cases  in  that  the  low  temperature 
chemical  rate  data  are  used. 

The  first  card  for  each  case  consists  of  a  set  of  integers  which  are 
indicators  used  in  tne  program.  This  card  must  always  be  first.  Thereafter, 
the  input  data  should1  follow  the  following  format  requirements.  There  are 
five'  data  entries  per  card  preceeded  by  the  location  of  the  first  item 
(Columns  1-9)  on  the  card.  This  format  facilitates  reading  a  variable  number 
of  pieces  of  real  data  into  a  specified  location  in  an  array;,  as  a  result  the 
order  of  these  data  cards  to  one. another  is  immaterial  and  they  may  be  in 
sequential  or  nonsequential  order. 

The  number  (index)  in  positions  1-9  of  the  first  field  on  each  card  de-  , 
fines  the  location  of  the  first  piece  of  data  on  that  card*  If  the  index  is 
101,  the  first  piece  of  data  is  stored  in  the  101st  location;  the  remaining 
fields  on  each  card  contain  data  for  the  successive  locations.  If  one  or  more 
of  the  data  fields  are  left  blank,  no  information  is  read  into  the-  locations 
corresponding  to  these  fields;  the  information  already  in  these  locations 
(if  any)  remains  unaltered.  Hence,  when  running  successive  cases,  certain 
items  may  be  changed  leaving  others  as  they  were  previously. 

Data  read  using  the  method  described  above  must  conform  to  the  follow¬ 
ing  specifications : 

1.  The  index  must  be  written  to  the  extreme  right  of  positions  l-9> 

and  this  index  must  not  be  zero  or  blank.  >  - 

2,  The  floating  point  (real)  data  written  in  the  five  data  fields 
may  be  written  in  two  ways: 

(a)  with  a  decimal  point  and  no  exponent 

(b)  with  the  E12; 8  type  format  where  the  decimal  point -is 
located  eight  places  to  the  left  of  the  beginning  of  the 
exponent  field, 

3 .  Data  reading  is  concluded  by  placing  a  negative  sign  ahead  of  the  index 

of  the  last  card  for  the  case.  Data  information  may  or  may  not  be  in 

the  fields  of  this  card. 
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Sines  the  input  region  is  not  cleared  between  cases,  the  input  data  from 
a  previous  case  roay  be  used  in  subsequent  cases,  thus,  when  several  cases 
are  to  be  run  consecutively  and  some  of  the  input  data  are  the  same,  these 
data  items  need  not  be  re-entered, 

5.  PROGRAM  OUTPUT  - 

A  complete  set  of  sample  printout  is  shown  in  Appendix  III  for  the  HAM- 
83  streamline  No.  5  case.  The  first  hatch  of  printing  id  the  input  data  for 
the  case.  Again,  all  FORTRAN  symbols  are  defined  at  the  beginning, of  the 
report  and  for  the  most  part  the  output  is  self-explanatory.  The  term. 

"entries.  1-5,"  etc.,  refers  to  the  order  in  which  the  chemical  reactions  „  „ 

appear  in  the  program.  The  maip  printout  presents  the  local  streamtube  co-  .  „„ 
ordinates  X  and  Y  and  the  streamtube  distance,  S,  as  measured  from  the  start¬ 
ing  point.  Next,  the  gross  flow  field  properties  (pressure,  velocity,  stream- 
tube  area,  density,  enthalpy,  and  translational  temperature)  are  tabulated. 

The  parameter  EETT  is  included  for  reference  only.  It  was  originally  the 
electronic  energy  based  on  translational  temperature.  It  has  recently  been 
changed  to  the  electronic  energy  based  on.  the  vibrational  temperature  of  nitro¬ 
gen.  Following  this  output  are  the  chemical  composition,  vibrational  temper¬ 
atures,  electronic  temperatures,,  the  vibrational  energies,  and  the -electronic 
energies . 

Between  the  printout  for  each  integration,  interval  one.  observes  spme 
miscellaneous  information.  This  information  is.  basically  for  diagnostic 
purposes  and  is  purposely  printed  on  separate  pages  so  that  it  can  readily  . 
be  detached  and  eliminated  from  the  main  printout  if  so  desired.  The  .  eras 
labeled  HARMON,  DES,  and  FORM  represent  the  time,  constant  term,,  the  destruc¬ 
tion  term,  and  the  formation*  term,  respectively.,  in  the  vibrational  rate 
equation.  The  other  printout  applies  to  the  stability  and  accuracy  tests  in 
the  integration  subprogram.  It  indicates  which  integrated  parameter  is  set¬ 
ting  the  integration  step  size  and  under  what  conditions  it  is  doing  so.  The 
line  of  fixed  point  numbers  represents  the  difference  between  the  ordinates 
Y2  and  Y3  (i.e.,  the  denominator  in  the  stability  test)  so  that  one  can  see 
which  parameters  would  have  denominators  lying  below  the  acceptable  value  of 
75.0.  This  applies  to  the  accuracy  of  the  stability  test. 

The  order  in  which  the  chemical  species  are  considered  both  in  the  pro¬ 
gram  computation  and  in  the  program  output  is  always  the  same  and  as  shown 
below: 


N2,  02,  NO,  N0+,  N2+,  02+,  02“,  N,  0,  N+,  0+,  0",  e,  Ar 
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SECTION  IV 

RESULTS  FOR  SAMPLE  CASES 


>  f 

'  'N  •>  V 

1.  NORMAL  SHOCK  FLOW 

<■  % 

Preliminary  check  runs  were  made  for  the  case  pf  pureo^ygen  dissociation 
behind  a  normal  shock  including,  finite  vibrational  excitation. ’  The  results 
are  compared  with  those  obtained  by  Treanor  and  Karrone,  Reference  (2),  in 
Figure  8  .  Two  cases  were  run  with  the  S  J)  program.  The  first  used  the 
Millikan  and  White  vibrational  relaxation  time  constants,  Reference  (12),  with 
electronic  energy  .included  in  the  species  enthalpies.  In  the  second  run,.,  , 
Plackmn's  vibrational  relaxation  time  constants,,  also  used  by  Treanor  and  „ 
Marrone,  were  used  and  no  electronic  energies  were  considered.  „Xn  all  other  l 
respects,  both  the  formulation  and  data  corresponded  to  the  Treanor  and 
Mdrrcne  analysis.  The  two  SD  runs  correctly  show  the  expected  difference 
using  the  shorter  time  constants  which  result  from  using  the  Millikan  and 
White  correlations.  The  second  run  was  expected  to  duplicate  the  Reference 
ti“)  results,  but  as  can  readily  be  seen  in  Figure  S  *  produced  slightly  dif¬ 
ferent  results.  Although  the  temperatures  in  the  '’quasiequilibrium"  zone 
downstream  of  S  =  0.1  cm  appear  to  be  in  fair  agreement*  a  significant 
difference  is  apparent  in  the  initial  rate  of  vibrational  temperature  increase 
ir.  the  region  S  <  0.05  cm.  To  verify  the  SD  results,  exact  hand  calcula¬ 
tions  were  made  of  the  initial  $  -derivatives  of  all  quantities  being  inte¬ 
grated.  and  the  S  -derivative  of  O2  vibrational  energy  at  S  =  0.0232  cm.  The 
ir.iti&l  derivative  check  showed  exact  agreement.  At  S  -  0.0232  cm  the 
following  vibrational  energy  rate  term  was,  computed  for  the  specie  O2. 


*  glk  — 

*  L  ^  J 


4.42  x  10^  ergs/gm-qm  Hand  Calculation 
4.42  x  1011  ergs/gm-cm  S&ID  Program 
8.10  x  10^-  ergs/gm-cm  Treanor' and  Marrohe 


d  S 


This  particular  term  is  by  far  the  largest  one  in  determining  the  net.  effect, 
the  chemical  rate  terms  being  small  at  this  comparatively  early  point  In  the 
relaxation  process  behind  the  shock.  Hand  calculations  were  also  made  of  the 
initial  derivatives  of  all  variables  immediately  behind  the  "shock.  Program 
results  also  checked  the  hand  calculations.  It  is  concluded  that  the  SD 
program  is  correct. 

The  program  was  next  used  to  compute  the  same  case  with  preferential 
dissociation  using  \J  =  ,  a  value  recommended  in  Reference  3  as  being 

in  agreement  with  experimental  results.  The  effect  of  preferential 
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Distance  Along  Streamtube,  s  ^  cm 

Figure  8*  Normal  Shoe!-  Relaxation  in  Pure  Og 
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Several  normal  shock  cases  for  pure  nitrogen  were  run  with. the  Bortner 
chemical  model  and  Mlllfken  and  White  vibrational  relaxation  data.  These  runs 
are  described  mere  fully  in  Section  II-6  .  The  runs  included  electronic 
excitation  effects.  It  is  of  interest  to  £hotf  the  data,  obtained  on  the  fina* 
such  run  illustrating  relative  Cw*ioribution  to  the  vibrational  energy 
of  »2  fi'on  the  relaxation  term  arid  the  chemical  formation;  and  destruction 
tern®.  This  comparison  is  presented  in  Figure  10. 
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2  HYPERSONIC  NOZZLE  FLOW 


A  hypersonic  nozzle  flow  case  simulating  conditions  in  the  three  megawatt 
facility  at  the  Air  Force  Flight  Dynamics  laboratory  (AFFDL)  was  run  with  the 
SD  program.  The  tunnel  reservoir  conditions  for  this  case  were  specified  as 

PT  =  132.68  ATM 

Ht  *  4284.76  BTU/LBm 

Tt  =  5725. 1°K 

The  starting  point  data  along  with  a  sample  printout  of  a  nozzle  run  made 
with  the  original  program  described  in  Reference  1  was  furnished  SD  by  AFFDL. 
The  starting  point  was  taken  slightly  downstream  of  the  nozzle  throat  section 
and  the  flow  was  assumed  to  be  in  both  thermal  and  chemical  equilibrium.  Vi¬ 
brational  and  electronic  temperatures  were  set  equal  to*  the  translational 
temperature. 

The  selection  of  applicable  chemical  species  for  the  air  model  is  based 
'on  the  species  that  would  be  present  ii  mole  fractions  greater  than  ICT^  for 
either  the  nozzle  case  dr  the  RAM-B3  case.  Reactions  14,  15,  and  20  were 
included  for  completeness  in  comparing  RAM-B3  results  even  though  the  ionized 
species  are  not  present  with  mole  fractions  greater  than  10~6.  These  are  the 
dominant  ionization  reactions  for  these  cases. 

Input  sheets  for  the  nozzle  case  are  presented  in  Appendix  XI.  Included 
with  the  initial  starting  point  data  is  a  description  of  the  area-match  para¬ 
meter,  TABAP,  as  a  function  Of  nozzle  station*  TABX.  This  latter  data  is 
based  qn  the  nozzle  contour  defined  in  Reference  22,  .Again,  it  should  be 
emphasized  that  the  input  data  presented  in  Appendix  XI  represent  changes  to 
the  sample  data  for  the  RAM-B3,  streamline  No.  5  case. 

Figure  11  shows  the  results  of  the  SD  run  and  gives  the  results  of  the 
AFFDL  rim  for  comparison.  In  the  interest  of  simplicity,  only  curves  for  N 
and  0  are  shown  in  Figure  11,  The  electronic  temperature  for  N2  appears  to 
freeze  immediately  whereas  a  delay  in  the  freezing  is  experienced  for  the 
vibrational  temperature  of  N2.  The  initial  mole  fractions  of  N  and  0  were, 
identical  for  the  two  runs  arid  the  decrease  in  the  mole  fraction  of  N,  shown 
in  Figure  11,  takes  place  immediately.  It  can  probably  be  traced  to  alight 
differences  in  the  chemical  rate  data  for  the  two  programs. 

Sample  output  sheets,  representative  of  the  results  shown  in  Figure  11, 
are  also  presented  in  Appendix  XI.  The  number  of  integration  intervals 
(ICOUNT)  and  hence  the  run  time  proved  quite  large  for  the  nozzle  case.  This 
is  partially  due  to  starting  the  problem  under  near-equilibrium  conditions 
where  the  reaction  equations  are  expected  tb  be  ’'stiff”  arid  integration  stab¬ 
ility  is  indeed  a  major  factor  in  setting  step  size.  More  important,  however, 
was  the  fact  that  the  program  was  still  in  checkout  during  the  running  of  the 
nozzle  case  and  five  separate  program  runs  were  required  to  complete  the 
nozzle  case.  Over  one  segment  (the  longest  run)  the  integration  procedure 
was  not  properly  increasing  step  size.  Ibis  deficiency  was  corrected  on  the 


subsequent  runs.  The  actual  run  took  thirty-five  minutes, Ip  complete?  but  it 
is  felt  that  with  the  final  version  of  the  computer  program  the  run  time  would 
only  require  about  fifteen  to  twenty  minutes.  It  should  be  noted  that  while 
the  program  for  this  run  was  inefficient,  the  results  are  considered  accurate. 

3.  RAM  53  FLOW  FIELD  ANALYSIS  '  ’ 

The  pressure  matching  portion  of  the  new  program  was  used  to  generate  the 
external  nonequilibrium  flow  field  for  the  RAM  B3  reentry  configuration.  The 
flow  field  generated  was  for  Case  No,  2  in  the  RAM  B3  report,  Reference  23  , 
where  the  flight  conditions  are  as  Jb'iiows:  ^ 

Velocity  =  10,.7Q7  fps  .  . .  ? 

Altitude  «  135,539  [ft 

The  RAM  B3  configuration  is  4  spherically-blunted,  nine-degree,  semi-apex 
argle  cone  which  is  shown  in  Figure  12  along  with  the  eight  streamlines  for 
which  data  were,  generated  during  the  course  of  this  contract..  The  streamlines 
selected  were  direct  counterparts  of  streamlines  analyzed  in  Reference  23. 
Figure  12  presents,  in  addition  to  the  selected  streamline  pattern,  a  series, 
of  body  normals  for  which  data  is.  presented  in  Reference  23  •  The  dashed  line 
shown  in  the  figure  represents  the  boundary  layer  build-up  over  the  body  sur¬ 
face. 


Results  were  obtained  for  eight  streamlines,  and  data  plots  are  presented 
ir-  Figures  13,  14  ,  and  15  for  three  typical  streamlines  (streamlines  No.  5, 
13,  and  26).  Beyond  streamline  No.  26  very  little  chemical,  vibrational*  arid 
electronic  energy  excitation  occurs.  All  three  figures  have  the  positions  Pf 
normals  designated  so  that  with  the  aid  of  Figure  12  one  can  readily  identify 
trends  in  the  results  with  position  along  the  streamline. 

Streamline  No.  5  was  selected  as  the  typical  input-output  case  in  Section 
III.  It  was  also  the  only  one  of  the  .sight  to  enter  the  boundary  layer  and 
all’ computation  was  stopped  at  that  point.  Changes  in  input  data  and  sample 
output  results  for  the  other  streamlines  are  furnished  in  Appendixes  IV  -  X. 
The  matching  parameter  history  and  the  coordinates  of  the  streamline  are  a 
part  of  the  input  4ata.  Th0  run  time  for  these  streamlines  varied  from  one 
minute  on  streamline  No.  38  to  about  five  minutes  on  streamline  No.  5* 

Starting  point  conditions  for  RAM  B3  streamlines  were  obtained  from  shock 
calculations  assuming  the  translational  and  rotational  energy  modes  are  im¬ 
mediately  equilibrated  behind  the  shock.  The  chemical,  vibrational,  and  elec¬ 
tronic  energs  modes  on  the  other  hand  are  assumed  to  reniain  at  free  stream 
conditions  through  the  shock  itself  and  then  proceed  to  relax  as  the  flow 
mores  downstream. 


The  initial  concentrations  of  species  N  and  0  were  set  at  10*^  rather 
than  10“°,  as  was  the  case  in  Reference  23,  to  give  a  more  realistic  represen¬ 
tation  of  the  iiiitial  transients  behind  the  shock.  The  results  for  streamline 
iji 0.  26  indicate  even  lower  initial  concentrations  of  these  species  are/ desirable 
but  the  program  requires  that  all  species  included  in  the  model  be  present  in 
finite  quantities  for  it  to  operate  properly. 
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Streamline  Pattern. 


Data  from  Ref,  23  (Tt  *  ?r  *»  Te) 
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Figure  14.  RAM-B3  ^equilibrium  Flow  Properties  s 
Streamline  No.  13 
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APPENDIX  I 


.  VIBRATIONAL  ENERGY  BALANCE  . 

•.  ■  ;  „,r  jY: 


■.  ?  -<• 


The  derivation  of  the  vibrational  energy  equation  (Equation  (8)  in  the 
main  text )  is  based  bn  a  simple  energy  balance.  Consider  d&vtfat  to  be  the 
time  rate  of  change  of  vibrational  energy  per  gram  of  mixture.  This  rate  is 
related  to  two  effects:  (l)  the  vibrational  energy  transfer  due  to,  molecular 
translation  effects  (no  chemical  effects),  and  (2^ the  effect  on  vibrational 
energy  resulting  from  chemical  changes.  This  can  be  put  into  equation  form 
as  follows:  , 


dtrvSl  -  1  time  rate  of  vibrational  energy  transfer,  due  to  translation 

TE 


-E1 

1 

f 


g  (mixture) 

vibrational  energy  lost\  /molesi  dissociated. 

1Jt  _  \  seC-g  (mixture),  J\ 

!  -  '  '  l..  -  -  <  .f 

+ j ^vibrational  .energy'  gained V molesj  recombined  \ 
molef  aec-g  (mixture)  / 


1 

m 


Substituting  the  physical'  -terms  Into  Equation  (i-l)  yields 

*  .  o  .  ”  **  '  ■» 

Mu  Nx 


cfcnr; 

\  M  / 


HtUt 


Me 


The  quantities  E»J£»o  »  >  and  1  are  evaluated  from,  the  energy  ex¬ 

pression  for  the  cutoff  harmonic  oscillator  model  shown  below  using  T:  ,  T5~c> 
and  Tf-j  ,  respectively.. 


The  temperature,  ,  is  defined  as 
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(1-3) 
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The  energy  term  GlT^.  is  obtained  from  Equation  (1-3)  using  T "V 

Since  we  are  basically  interested' in  the  vibrational  energy  rate  , per 
gram  of  species  £  s  for  use  in  Equation  (18),  we  can  transform 
(1-2)  by  means  of  the  relationship 


6<f.  - 


d-5)  - 


Differentiating  Equation  (1-5)  yields 


fori  *  —  Jmc  4-  t  d  iSVy 
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(1-6) 


Substituting  (1-2)  into  (1-6)  and  changing  from  time  to  streamline  distance 
as  the  independent  variable  produces  the  relationship 


which  is  identical  to  Equation  (8)  in  the  main  text 


APPENDIX  II  i 

VIBRATIONAL  RELAXATION  IN  A  MULTIOOHPQ^NT  MIXTURE 


The  vibrational  relaxation  for  a  diatomic  speciee  Which  has  translational , 
collisions  with  a  series  of  catalytic  specj.es  arid  i3  vibrationally  excited  by  *• 
eacn  encounter  must  account  for  each  separate  collision  in  its  overall  makeup.. 
Consider  a  single  molecule  of  species  A  and  the  possibility  that  it  will  en¬ 
counter  either  other  A*  molecules  or  molecule*.  of  species  B  or  C.  The  vibra¬ 
tional  energy  exchanged  in  the  collision  process.  is  obtained,.  in,  the  following 
manner: 
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Each  individual  vibrational  energy  rate  term  in  Equation  £±T*±).  is  .  defined  as 
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It  should  be  noted  in  this  analysis  that  the  first  of  the  two  reacting  species 
(i.e.,  species  "A")  is  considered  present  in  high  dilution  in  the  second-named 
gas  and  for  this  reason  each  value  of  is  based  on  the  partial  pressure 

of  the  second  gas  media.  The  energies  in  the  equations  presented -above  are 
defined  in  terms  of  energy  per  molecule /of  species  A.  Substituting  Equation 
(II-2)  into  Equation  (±1-1)  yields 
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where  'f*  represents  the  mixture  of  gases  and  is  defined  as 
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, _  f  t  it  it 
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(11-42 


The  expression  for  the  individual  values  of  tr  are  obtained  from  the  .relation¬ 
ship  -  -  '•  ,  V,  /:  <  . 
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um*-*  (Ti-e.oiff^  |-fe.4*|  ,  (n^) 


which,  when  substituted  in  Equation  yields 
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The  relationship  given  in  Equation  (II— 3)  can  be  generalized  to  apply 
to  the  complete  gas  media  by  multiplying  each  term  in  the  equation  by  the 
number  of  molecules  of  species  A  per  gram  of  mixture  producing 
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Each*£  is  still  baaed  on  the  partial  pressure  of  the  "J"  constituent. 
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In  iL\v  mixture  of. perfect  gases,  Dalton's  Law  states  that  ^ 
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It  is  easily  shown  that 


wh^re  denotes  mss  and  X  represents  the  mole  fracfciqn  of  the  specif  In 
qiestion^^Redefiidng  't  for  each  constituent  in  terms  of  the  total  j^essure, 
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‘ Combining  Equation  (11-6)  4nd  (11-9)  one  obtains  1 
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and  applying  the  mole  fraction  criteria,  of  Equation  (II-8)  one  finally 
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where  ^  ^  ,  etc,  ,  are  all  based  .on. the  pressure  of  themaix- 

ture.  IMs  is  the  result  derived  in  Reference  12  and  shown  in  Equation,  (30)* 
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<>■  ABSTRACT 


A  computer  program  has  been,  developed  for  one-dimansional  nonequillbrium 
reacting  gas  flow.  The  program  is  written  iii  Fortran  IV  and  is  compatible  with 
the.  IBM  7044/7094  direct  coupled  digital  computer  system  at  Wrighfc-Pattarfion  Air 
Force  Base,  Ohio.  In  addition  to  nonequilibriua  chemistry,  the  program  includes 
nonequilibrium  vibrational  and  electronic  energy  xelaxationand  coupling  effects 
between  these  energy  modes  and  the  chemistry.  The  formulation  is  based  on  a  one- 
dimensional  flow  matching  either  a  prescribed  pressure  car  area  variation  along  a 
streambube.  Therapdynamic  properties  are  computed  by  assuming  an  ideal.,, gas  mix¬ 
ture  and  the  equilibration  of  translational  and  rotational  temperatures.  ‘Rue 
jj  internal  energy  modes,,  rotation,  vibration,  and  electronic  excitation,  are  con¬ 
sidered  uncoupled;  and  a  rigid  rotator,  cut  off  simple  harmonic  oscillator, 
independent  of  the  electronic  state,  is  assumed.  Excitation  of  vibrational  and 
electronic  energies  are  treated  similarly  with  terms  which  account  for  relaxa¬ 
tion  and  chemical  reactions.  The  effects  of  nonequilibrium  vibrational  and 
electronic  states  on  chemical  rates  are.  included  ip  the  coupling  analysis.  The  . 
vibrational  relaxation  time  constants  were  obtained  from  the  KilHfcan  and  White 
data  while  the  electronic  relaxation  -tins-  constants  were  determined  for  nitrogen 
from  an  analysis  of  existing  shock  tube  radiation  measurements.  The,  computer 
program,  was  used  to  solve  for  the  nonequilibrium  flow  in  a  hypersonic  nozzle  and 
for  eight  streamlines  in  the  invlscid  flow  field  over  a  spherically  blunted  nine- 
degree  semiapex  angle  cone  at,  zero  angle  qf  attack. 
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